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ABSTRACT 


A  study  has  been  made  to  relate  theoretically  predicted  and  experimentally 
observed  pressure-coupled  combustion  instability  behavior  of  solid  propellants. 
Using  a  generalized  theoretical  model  which  incorporates  both  surface  and  gas 
phase  combustion  effects,  the  linear  and  nonlinear  (higher  pressure  amplitude) 
region  has  been  characterized  by  a  group  of  combustion  parameters.  Concurrent 
experimental  studies  have  been  directed  toward  investigating  (1)  the  effects 
of  coating  the  solid  oxidizer  crystals,  (2)  the  chemical  and  physical  nature 
of  the  propellant  binder,  catalysts,  and  cross-linking  agents  on  the  acoustic 
admittance,  (3)  nonlinear  response  properties,  and  (4)  propellant  burning  rate. 
The  data  show  that  small  changes  in  propellant  formulations  can  have  a  signifi¬ 
cant  effect  on  the  derived  parameters  characterizing  the  surface  and,  therefore, 
on  relative  contributions  of  these  reactions  to  the  acoustical  response  proper¬ 
ties  . 


A  comparison  of  the  experimental  data  with  theoretical  predictions  in  the 
linear  region  (i.e.,  acoustic  response  functions)  shows  that  propellants  having 
84%  solids  loading  are  more  likely  to  show  the  surface  gasification  to  be  pres¬ 
sure  dependent  than  those  having  lower  loadings.  The  temperature  dependence 
of  the  gasification  process  is  extremely  large —  larger  than  can  be  explained 
by  an  Arrhenius  rate  expression  with  a  reasonable  activation  energy.  The  rate 
of  energy  transfer  into  the  solid  phase  was  found  to  increase  with  increasing 
surface  temperature.  This  suggests  the  influence  of  exothermic  surface  processes. 
The  effect  of  pressure  on  the  energy  transfer  was  found  to  be  consistent  with  the 
severa\  proposed  models  for  gas-phase  combustion  effects. 

Higher  order  effects  show  that,  depending  on  the  magnitude  of  the  constants, 
both  positive  and  negative  values  for  the  shift  in  the  time  average  burning  rate 
with  acoustic  pressure  can  be  obtained.  Since  Eisel  has  observed  this  effect 
to  be  negative,  this  defines  the  region  of  reasonable  second  order  coefficients. 
Additional  studies  show  that  the  acoustic  response  function  decreases  with  in¬ 
creasing  acoustic  pressure  for  reasonable  second  and  third  order  coefficients. 
Hence,  the  combustion  process  itself  can  limit  the  acoustic  pressure  amplitude 
in  the  T-burner.  On  the  basis  of  these  calculations,  limiting  amplitudes  of 
0.1  to  0.5  can  be  expected  without  consideration  of  acoustic  losses. 

The  combined  linear  and  nonlinear  results  indicate  that  a  homogeneous  non¬ 
reacting  solid  characterizes  the  transient  behavior  of  the  solid  phase  with 
sufficient  accuracy.  Comparison  of  the  predicted  response  function-frequency 
curve  with  the  experimental  data  shows  agreement  within  reasonable  experimental 
errors.  The  nonlinear  behavior  also  shows  that  predicted  results  are  consistent 
with  observations.  Although  more  complex  and  detailed  descriptions  can  be 
proposed,  this  assumption  represents  an  extremely  useful  and  sufficiently  accu¬ 
rate  approximation  for  most  practical  purposes. 
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1.0  INTRODUCTION 

Unstable  combustion  of  solid  rocket  propellants  is  of  interest  for  two 
reasons:  (1)  it  is  one  of  several  methods  available  for  investigating  the 

structure  of  the  solid  propellant  combustion,  and  (2)  unstable  combustion 
has  presented  serious  problems  in  the  development  of  some  operational  solid 
propellant  rocket  motors.  Under  Contract  No.  DA-04-200-AMC-968(x) ,  UTC 
has  been  investigating  certain  aspects  of  this  problem,  namely  pressure- 
coupled  combustion  instability.  The  principal  objective  of  these  studies 
has  been  to  examine  the  relationship  belween  combustion  instability,  as 
observed  by  pressure  oscillations  in  thi  gas  phase,  and  exothermic  reactions 
occurring  on  and  within  the  solid  phase  surface. 

Experimental  studies  have  been  directed  toward  investigating  (1)  the 

effects  of  coating  the  solid  oxidizer  crystals,  the  chemical  and 
0 

physical  nature  of  the  propellant  binder,  catalysts,  and  cross-linking 
agents  on  the  acoustic  admittance,  (2)  nonlinear  response  properties,  and 
(3)  propellant  burning  rate.  Concurrent  theoretical  studies  have  been 
directed  toward  predicting  the  effects  of  these  variables  on  the  acoustic 
response  function.  Additional  analytical  studies  were  conducted  to  study 
the  transition  from  linear  to  nonlinear  acoustic  response  as  the  amplitude 
of  the  pressure  oscillations  increase.  This  report  represents  a  summary 
of  the  work  accomplished  under  Contract  No.  DA-04-200-AMC-968(x) . 
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2.0  TECHNICAL  DISCUSSION 

Recent  studies  on  the  structure  of  the  combustion  zone  of  composite 

solid  propellants  indicate  that  there  are  exothermic  processes  which  occur 

Cl  through  7)* 

on  and  within  the  solid  phase.  The  results  of  these  studies 

indicate  that  the  combustion  process  is  controlled  by  two  interdependent  exo¬ 
thermic  reaction  zones  near  and  on  the  surface  of  the  propellant.  One  zone  is 
in  the  gas  phase  at  a  finite  distance  from  the  solid  propellant  surface  and  is 
characterized  by  interdiffusion  of  gasified  oxidizer  and  fuel  species  and 
combustion  of  particles  of  ejcctrd  matter  from  the  surface.  The  second  reaction 
zone  occurs  on  and  within  the  solid  propellant  surface.  The  primary  heat  re¬ 
lease  in  this  zone  can  occur  as  a  result  of  exothermic  decomposition  processes 
and  chemical  reactions  between  the  initial  decomposition  products  of  the  various 
propellant  constituents.  Transient  and  steady-state  combustion  studies  indicate 
that  part  of  the  pressure-dependent  combustion  process  can  be  associated  with 
these  reactions. 

The  exothermic  surface  reactions  can  release  sufficient  heat  to  expel 
partially  combusted  products,  pyrolysis  products,  and  fuel  and  oxidizer  frag¬ 
ments  into  the  gas  phase  zone  above  the  surface,  where  they  intermix  and  burn 
completely.  The  maximum  flame  temperature  is  reached  in  the  luminous  zone, 
where  the  largest  portion  of  the  heat  is  released.  However,  because  of  the 


Parenthetical  superscript  numbers  denote  references  appearing  on  page  80. 
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relatively  large  mass  flow  perpendicular  to  the  surface,  only  a  small  amount  of 
heat  released  in  the  luminous  flame  zone  can  normally  reach  the  surface  to  sup¬ 
plement  the  heat  generated  by  the  surface  and  subsurface  reactions. 

Incorporation  of  the  exothermic  chemical  processes  at  the  solid  surface 
represents  an  important  addition  to  the  analysis  of  propellant  combustion 

/O  UN 

phenomena.  Until  recently,  theoretical  treatments  of  steady-state  combustionv  ' 
as  well  as  combustion  instability v  have  considered  the  effect  of  the  exo¬ 

thermic  combustion  reactions  to  originate  predominately  in  the  gas  phase.  How¬ 
ever,  as  Friedly^^  and  Marxman^"^  have  shown  theoretically,  the  principal 
time  lags  in  the  combustion  zone  are  associated  with  energy  transport  within 
the  solid  phase  for  oscillations  below  10,000  cps.  By  comparison,  the  time 
constants  associated  with  the  reaction  and  transport  processes  in  the  gas  phase 
are  small,  and  these  processes  can  be  considered  to  be  in  equilibrium  at  any 
instant.  Because  the  reactions  on  and  below  the  surface  are  exponentially 
dependent  on  the  solid  phase  temperature,  the  energv  released  by  these  reactions 
should  be  an  important  consideration  in  the  analysis  of  transient  combustion 
phenomena.  Furthermore,  the  reported  data  suggest  that  these  reactions  exhibit 
a  significant  dependence  on  pressure.  Thus,  the  acoustic  pressure  oscillations 
will  produce  oscillations  in  reaction  rates  through  both  the  pressure  and  tempera¬ 
ture  sensitivities  of  the  reaction  rates  which  can  have  a  significant  effect  on 
the  combustion  stability  behavior  of  solid  propellants. 

2 . 1  EXPERIMENTAL  STUDIES 

The  experimental  studies  have  been  directed  toward  investigating  the 
effects  of  propellant  composition  variables  on  the  response  of  the  combustion 
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process  to  acoustic  pressure  oscillations.  The  principal  propellant  property 
derived  from  these  studies  is  the  real  pare  of  the  acoustic  response  function, 
Re(u/e),  which  relates  the  amplitude  of  the  burning  rate  oscillations  to  the 
amplitude  of  the  incident  acoustic  pressure  oscillations  by  the  relationship 

(r/h  ) 

Re(y.)  =  — —  (l) 

e  <p/po) 

where 

r  =  steady-state  turning  rate 

P  =  mean  combustion  pressure 

o 

r  =  amplitude  of  the  burning  rate  oscillations 
¥  =  amplitude  of  the  acoustic  pressure  oscillation 

The  approach  taken  in  these  studies  has  been  to  change  the  propellant 
composition  in  such  a  way  that  contributions  from  the  possible  exothermic  sur¬ 
face  reactions  are  varied  and  to  observe  the  resultant  effects  on  the  acoustic 
response  function. 

An  attempt  was  also  made  to  examine  the  propellant  combustion  visually 
under  oscillatory  conditions.  A  series  of  preliminary  experiments  was  conducted 
in  which  the  combustion  process  was  altered  by  phase  shifting  the  pressure 
oscillations,  These  topics  are  discussed  in  appendix  A. 

2.1.1  Experimental  Apparatus 

Measurement  of  the  acoustic  response  function  can  be  obtained  in  a  T-burner 
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apparatus.  This  type  of  equipment,  which  has  become  almost  standard  in  the 

industry  for  acoustic  measurements,  consists  of  a  cylindrical  combustion  bomb 

with  end  burning  charges  of  propellant  in  cither  or  both  ends.  By  measuring 

the  oscillating  pressure  at  the  ends  of  the  burner,  it  is  possible  to  derive 

the  acoustic  response  furction  fsr  the  particular  propellant  under  study.  The 

(19) 

apparatus  used  in  these  studies  is  based  on  the  design  concepts  of  Horton 
and  Strittmater^^and  is  shown  schematically  in  figure  1.  A  photograph  of 
the  equipment  appears  in  figure  2. 


81309 


Figure  1.  Diagram  of  T-Burner 
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The  combustion  chamber  has  an  internal  diameter  of  1.5  in.  over  the  entire 
length  and  can  be  operated  as  either  a  single-  or  double-ended  T-burner.  Pro¬ 
pellant  samples  are  located  in  the  end  caps,  which  also  contain  provision  for 
the  igniter  leads  and  pressure  transducers.  Nichrome  ignition  wires  are  attached 
to  two  electrodes  which,  in  turn,  are  located  in  a  ceramic  tube  inside  the  cap 
and  pressure  sealed  with  a  connax  fitting.  The  nichrome  wire  is  passed  through 
the  propellant  sample  and  is  secured  to  the  surface  with  small  staples.  Powdered 
propellant  is  then  placed  over  the  wires.  Ignition  is  obtained  by  electrically 
heating  the  nichrome  wire  secured  to  the  propellant  surfaces.  The  two  end  caps 
are  shown  with  and  without  a  propellant  sample  in  figure  3.  The  end  cap  with  the 
propellant  sample  also  shows  the  ignition  wires.  Combustion  products  and  pres¬ 
surized  gas  are  exhausted  into  two  blowdown  tanks  through  four  orifices  equally 
spaced  around  the  circumference.  The  disturbances  and  resulting  acoustic  losses 
caused  by  this  secondary  flow  are  minimized  by  locating  the  orifices  at  the 
pressure  node. 

Both  the  combustion  tube  and  blowdown  tanks  can  be  prepressurized  with  any 
inert  gas  to  any  pressure  from  0  to  1,000  psig.  A  bleed  valve  is  used  to  exhaust 
the  system  to  atmospheric  pressure,  and  the  excess  gases  which  remain  in  the 
system  are  evacuated  with  a  vacuum  pump.  The  frequency  of  the  acoustic  pressure 
oscillations  can  be  varied  by  using  interchangeable  sections  of  different 
lengths,  as  shown  in  figure  2.  Using  the  double-ended  configuration,  frequencies 
between  150  and  5,000  cps  can  be  studied  with  the  present  design.  This  represents 
the  region  of  interest  since  the  response  function  curves  appear  to  go  through 
a  maximum. 
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The  steady-state  and  oscillating  pressures  were  measured  by  Kistler  model 
No.  601  pressure  transducers  in  contact  with  the  back  surface  of  each  propellant 
sample.  The  charge  amplifiers  were  mounted  in  the  test  cell  next  to  the  trans¬ 
ducers  to  reduce  the  line  losses.  Both  the  average  chamber  pressure  and  the 
oscillating  component  of  pressure  were  recorded.  To  record  the  average  chamber 
pressure,  the  transducer  signal  passed  through  a  60-cps  filter  and  was  recorded 
on  a  CEC  recording  oscillograph  using  a  No.  7-346  galvanometer.  The  high- 
frequency  signal  from  the  transducers,  which  was  used  to  determine  pressure 
amplitude  and  wave  forms,  was  recorded  on  an  Ampex  tape  recorder.  The  tape  out¬ 
put  was  then  passed  through  a  Spenser-Kennedy  band-pass  variable  electronic 
filter  and  amplified  with  a  McIntosh  amplifier  which  provides  the  gain  and  power- 
necessary  to  drive  a  CEC  No.  7-361  high-frequency  galvanometer.  Zener  diodes 
were  used  as  voltage  cutoff  to  protect  the  galvanometer.  To  eliminate  any  low- 
frequency  signals  (below  150  cps)  from  influencing  the  oscillatory  pressure  trace, 
a  capacitor  was  installed  between  the  McIntosh  amplifier  and  the  galvanometer. 
Thus,  a  low-frequency  noise  was  removed  from  pressure  traces,  thereby  improving 
the  accuracy  of  the  measurements. 

A  schematic  of  the  electronics  is  shown  in  figure  4,  and  a  photograph  of 
the  main  components  a' pears  in  figure  5.  The  response  of  the  instrumentation 
was  investigated  with  an  oscillator  over  the  frequency  range  from  150  to  5,000 
cps.  The  recorded  wave  forms  and  frequencies  were  not  altered  by  the  signal 
processing  electronics. 
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Figure  4.  Schematic  of  T-Burner  Electronics 


At  low  frequencies  the  combustion  chamber  becomes  long,  and  the  quantity 
of  room  temperature  gas  required  to  pressurize  the  chamber  becomes  significant. 
After  ignition,  a  significant  period  of  time  is  required  to  achieve  thermal 
equilibrium  in  the  combustor.  The  acoustic  pressure  oscillations  tend  to  grow 
during  the  period  in  which  thermal  equilibrium  is  being  established.  This 
contrasts  with  the  situation  in  which  the  instability  propagates  from  a  steady- 
state  condition,  which  usually  occurs  in  actual  motors  and  which  is  treated 
theoretically.  To  obtain  more  realistic  measurements  in  the  low-frequency  range, 
four  pieces  of  propellant  1/4  in.  wide  by  3/4  in.  long  by  1/16  in.  thick  were 
mounted  in  the  center  of  the  tube  near  the  four  vent  ports.  Because  these 
propellant  pieces  were  located  at  a  pressure  node,  the  oscillations  tended  to 


10 


Figure  5.  T-Burner  Electronics 


UTC  2136-FR 


be  suppressed  until  these  small  samples  burned  out,  allowing  the  propellant 
samples  in  the  T-burner  end  caps  to  achieve  steady  state  before  the  oscillations 
grew. 

2.1.2  Propellant  Formulations 

There  are  several  potential  methods  for  controlling  the  contributions  of 
surface  reactions  to  the  acoustical  response.  One  method  involves  coating  the 
oxidizer  crystals  to  alter  the  reactivity  of  the  crystal-binder  interface.  Also, 
changes  in  the  chemical  structure  of  the  binder,  in  the  agents  used  to  crosslink 
the  binder,  in  catalyst  type  and  content,  and  in  oxidizer  particle  size  and 
loading  were  all  postulated  to  influence  the  acoustic  response  function.  Accord¬ 
ingly,  propellants  were  prepared  which  permitted  investigation  of  these  variables. 

The  first  phase  of  the  experimental  program  examined  the  effect  of  coating 
the  oxidizer  with  different  polymeric  materials.  A  composite  propellant  normally 
contains  some  fine-grind  oxidizer;  however,  coating  fine  particles  results  in 
particle  aggregation.  With  each  polymer  coating,  the  size  of  the  aggregated 
particles  vary,  and  the  time  involved  to  eliminate  this  problem  from  the  program 
proved  to  be  impracticable.  The  problem  does  not  occur  with  larger  particles; 
therefore,  AP  oxidizer,  having  Average  particle  size  of  190|i,  was  used. 

The  solve .it-nonsolvent  technique  was  used  to  apply  the  polymer  coat  to  the 
oxidizer,  'lhe  selected  polymer  was  dissolved  in  a  suitable  solvent,  and  the 
required  amount  of  AP  was  added  to  this  solution  in  a  mixer.  The  solvent  was 
then  evaporated  from  the  mix  until  the  gel  state  of  the  polymer  was  present. 

A  second  liquid  that  is  nonsolvent  for  the  polymer  was  then  added  to  the  polymer 
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solution  at  a  low  rate.  This  forced  the  polymer  out  of  solution  and  around  the 
particles  of  AP.  Enough  of  the  nonsolvent  was  added  to  harden  the  polymer. 
After  stirring  for  a  specified  length  of  time,  the  liquid  was  decanted  and  the 
coated  material  dried.  Table  I  identifies  the  polymers  and  additives  used. 


TABLE  I 
ADDITIVES 

1.  LiF  -  Lithium  fluoride 

2.  Cab-O-Sil  -  Silicon  dioxide 

3.  Kel-F-800  -  Polychlorotrif luoroethylene 

4.  Hypalon-30  -  Chlorosulfonated  polyethylene 

5.  Ethyl  cellulose  -  OH  groups  in  cellulose  partly  or 

completely  replaced  by  ethoxyl 
groups,  0~C2H5 

Samples  were  taken  from  each  batch  and  tested  for  a  coating  by  dispersing 
the  oxidizer  in  a  drop  of  water  on  a  microscope  slide.  The  slide  was  observed 
under  a  microscope,  and  as  the  AP  dissolved,  a  coating  was  observed  as  an  empty 
polymer  shell. 

A  series  of  propellants  was  prepared  using  the  coated  oxidizer  to  study 
the  effects  of  the  coatings  on  admittance.  The  use  of  190(a  AP  limited  the 
propellant  formulation  to  78%  oxidizer,  and  the  binder  system  was  a  carboxy- 
tetminated  polybutadiene  to  give  good  flow  properries  to  the  mix.  The  mixer 
used  for  preparing  the  formulations  was  a  1-pint  vertical  Baker  Perkins.  After 
mixing,  the  propellant  was  cast  into  special  molds  for  the  T-burner,  and  strands 
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were  made  to  acquire  the  burning  rate  of  the  formulation.  The  propellant  was 
cured  for  5  days  at  160°F,  The  formulation  containing  the  ethylcellulose-coated 
oxidizer  did  not  wet  properly  during  the  mixing  cycle,  but  the  propellant  cured 
in  a  normal  manner  after  the  propellant  was  placed  in  the  molds. 

Two  formulations  were  prepared  with  the  coating  material  added  to  the 
propellant  mix.  These  propellants  were  used  to  check  the  effect  of  the  location 
of  the  coating  (e.g.,  in  the  propellant  or  as  a  coating  on  the  oxidizer),  to 
evaluate  the  possibility  that  the  coating  could  be  stripped  from  the  oxidizer 
during  mixing,  and  the  effects  of  coating  porosity. 

The  particular  propellant  formulations  tested  and  the  resultant  burning 
rate  data  are  shown  in  table  II.  Burning  rate  measurements  were  made  in  a 
standard  strand  burner  at  three  pressure  levels:  200,  500,  and  1,000  psig. 

The  data  for  each  system  were  checked  at  the  200-psig  level  witn  the  total 
burning  time  measurements  in  the  T-burner,  i.e.,  burn  time  of  the  1/4-in. -thick 
samples.  The  difference  between  the  two  values  was  less  than  10%  in  all  cases. 
Because  of  the  experimental  scatter  in  the  burning  rate  measurements,  the  strand 
data  are  valid  to  within  ±10%  of  the  reported  values.  In  some  cases,  the  reported 
burning  rate  pressure  exponents  are  abnormally  high  at  200  psi,  especially  con¬ 
sidering  the  large  oxidiz-.r  particle  size  which  was  used.  However,  the  data  at 
higher  pressures  show  that  the  exponents  for  many  of  these  propellants  are  not 
constant  over  the  entire  range. 

It  should  be  noted  that  propellants  No,  16  through  No.  23  were  tested  for 
P  and  L*  extinguishment  behavior  by  Jensen/21^  The  catalysts  had  little  effect 
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on  the  L*  extinguishment  characteristics.  The  iron  oxide  also  had  no  influence 

on  the  P  behavior,  while  the  copper  chromite  increased  the  f*  requirement  by 

approximately  30%  for  the  0,5%  catalyzed  propellant.  In  addition,  the  PEP 

propellant  had  higher  P  and  L*  extinguishment  limits  at  a  fixed  pressure  in 

comparison  to  the  UTREZ  propellant.  The  use  of  larger  particle  AP  was  also 

found  to  increase  the  P  required  to  terminate  combustion  at  pressures  above 

0 

250  psi  and  to  decrease  the  P  required  below  250  psi.  At  200  psi,  the  larger 
particle  propellant  (400^)  has  a  critical  P  approximately  half  that  of  the  200^ 
propellant. 


2.1.3  Reduction  of  Experimental  Data 

The  relationship  between  the  observed  acoustic  pressure  oscillations  and 

the  acoustic  response  function  of  the  propellant  has  received  considerable 
(20  22  23  24) 

attention.  ’  ’  ’  Although  several  methods  have  been  developed  to  derive 

the  acoustic  response  function  from  acoustic  pressure  measurements,  the  only 
method  used  in  this  study  was  based  on  growch  and  decay  constant  determination, 
that  is. 


Re  (J) 


P 


o 

4cp  r 

o 


(2) 


During  the  experimental  program,  special  care  was  required  when  tests  were 
conducted  at  low  frequencies.  Under  these  conditions,  the  combustion  chamber 
became  long  and  the  quantity  of  room  temperature  gas  required  to  pressurize  the 
chamber  thereby  became  significant.  After  ignition,  a  significant  period  of 
time  was  required  to  achieve  thermal  equilibrium  in  the  combustor.  However, 
the  frequency  of  the  pressure  oscillations  depended  on  the  temperature  through 
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the  sonic  velocity,  i.e.,  c  ~  VT,  through  the  expression 

f  =  c/2L  (3) 

S 

The  sonic  velocity  in  equation  2  results  from  considerations  at  the  combus¬ 
tion  zone  and,  therefore,  should  be  evaluated  from  the  experimental  flame 
temperature.  The  approach  used  in  this  study  was  based  on  the  assumption  hat, 
at  high  frequencies  (short  combustion  chamber)  the  experimental  flame  tempera¬ 
ture  is  established  throughout  the  combus  tion  chamber  very  rapidly.  Thus,  c 

was  evaluated  from  equation  3  extrapolating  the  parameter  2Lf  to  high  frequencies. 

S 

This  amounts  to  a  partial  correction  of  the  generalized  acoustic  response  analysis 
(25) 

by  NOTS'  in  that  the  common  uncorrected  response  function  is  based  on  an 

average  speed  of  sound  (c  )  in  the  denominator  of  equation  2  instead  of  the  sonic 

o 

velocity  of  the  gas  in  the  combustion  zone  (c) .  This  report  suggests  that  the 

correction  factor  applied  to  the  data  should  be  based  on  a  correction  term  given 
2 

by  (c^/c)  .  Thus,  if  our  data  are  to  be  corrected  further  in  accordance  with 
their  results,  then  the  numerical  values  of  Re(p/p)  at  lower  frequencies  should 
be  multiplied  by  cq/c  and  would  be  reduced  by  a  small  amount  at  low  frequencies. 

The  frequency  at  which  the  maximum  occurs  would  not  be  altered  by  the  correction 
term. 


2.1.4  Experimental  Results 

Acoustic  response  function  of  propellants  shown  in  table  II  were  determined 

as  a  function  of  frequency  in  the  T-burner.  These  data  are  presented  in  figures 

2 

6  through  19,  using  the  normalized  frequency  8rtf^/rQ  as  the  parameter.  This 

(12)  (17) 

form  was  used  based  on  the  theoretical  studies  of  Hart  and  Friedly  and 
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gure  9.  Acoustic  Response  of  CTPB  and  UTREZ  Binder 
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Effect  of  Pressure  on  Acoustic  Response 


FREQUENCY 


UTX  -  9368:  84%  AP/16%  UTP.EZ  BINDER 

UTX  -  9715:  84%  AP  WITH  0.25%  Fe-O, 


Effect  of  Various  Coatings 
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the  theoretical  studies  to  he  discussed  in  the  following  sections.  These 

theoretical  studies  suggest  that  the  proper  form  ior  the  dimensionless  frequency 

should  include  the  propellant  thermal  diffusivity.  For  the  propellants  listed, 

-4 

this  propellant  property  should  be  relatively  constant  (approximately  3  x  10 
2 

in.  /sec);  therefore  it  has  been  omitted.  In  addition,  the  definition  of  the 
acoustic  response  function  is  such  that  it  should  approach  the  normal  burning 
rate  pressure  exponent  as  the  frequency  approaches  zero;  therefore,  the  data  are 
reported  showing  the  extrapolation  at  low  frequencies  to  the  burning  rate  pres¬ 
sure  exponent  obtained  from  the  burning  rate  measurements. 

For  completeness,  all  the  experimental  data  obtained  for  all  the  propellants 

shown  in  table  II  are  presented  in  appendix  B.  The  values  of  the  growth  frequency 

(f  ),  growth  rate  (a  ),  decay  frequency  (f ,)  ,  decay  rate  (a.)  and  Re(u/c)  are 
8  8 

reported. 


It  is  important  to  note  that  the  data  are  subject  to  significant  uncertainty. 

The  steady-state  burning  rates,  as  previously  indicated,  may  have  errors  as  high 

as  ±10%.  These  values  enter  into  the  calculations  of  the  response  functions, 

2 

i.e.,  p  ~l/ro>  and  the  frequency  functions,  8jrf^/rQ  .  In  addition  to  these 
uncertainties,  the  measurements  of  the  growth  and  decay  rates  are  subject  to 
an  uncertainty  of  ±10%;  therefore,  the  absolute  value  of  the  response  function 
u/e  is  valid  to  within  ±20%..  The  experimental  results  must  be  interpreted 
according ly . 

As  part  of  activity  associated  with  the  ICRPG  round  robin  group  on  combus¬ 
tion  instability,  the  acoustic  response  function  for  A-13  propellant  (an  NWC 
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formulation)  was  measured  at  UTC.  The  results  of  these  tests  agree  with  the 

(19) 

earlier  data  obtained  at  NWC  by  Horton,  as  shown  in  figure  10.  Hence,  it  is 
reasonable  to  assume  that  the  differences  in  the  design  of  the  two  T-burners 
do  not  influence  the  experimental  results. 

2.2  THEORETICAL  STUDIES 

Before  the  experimental  data  presented  in  the  preceding  paragraphs  can  be 
analyzed,  it  is  necessary  to  have  some  theoretical  framework  as  a  guide.  To 
develop  this  framework,  one  can  take  two  approaches.  For  the  first  approach,  a 
rather  complete  and  detailed  model  of  the  combustion  process  must  be  postulated. 
The  various  important  processes  are  characterized  mathematically  and  the  result¬ 
ing  equations  are  solved  to  obtain  the  desired  theoretical  result.  By  comparing 
theoretical  and  experimental  results,  one  can  derive  values  for  specific  physical 
constants  (i.e.,  activation  energies,  heats  of  reaction,  etc.)  which  are  presumed 
to  characterize  the  combustion  process.  The  results  reflect  the  particular  bias 
of  the  specific  model,  and  the  derived  values  for  the  various  physical  constants 
are  often  difficult  to  interpret  since  independent  measurements  are  difficult  to 
obtain. 

Another  approach  minimizes  the  specific  assumptions  concerning  the  detailed 
structure  at  the  combustion  zone.  The  result  is  a  mathematical  formulation 
which  contains  a  number  of  parameters  which  have  no  definition  in  physical  terms. 
They  can  be  derived  from  the  experimental  data  and,  by  examination  of  a  number  of 
possible  physical  definitions,  some  further  insight  into  the  combustion  zone  can 
be  obtained.  The  principal  advantage  of  this  approach  is  that  one  is  not  com¬ 
mitted  to  a  specific  model;  instead,  the  experimental  evidence  guides  the 
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theoretical  development.  Also,  by  limiting  the  number  of  physical  assumptions, 
one  can  better  test  their  validity. 

In  the  paragraphs  which  follow,  the  general  analysis  will  be  presented,  and 
then  the  specific  analysis  developed  early  in  the  program  will  be  derived  from 
the  general  treatment. 


2.2.1  General  Combustion  Analysis 

For  the  purpose  of  analysis,  the  solid  phase  within  the  combustion  zone  is 
assumed  to  be  characterized  by  the  transient  one-dimensional  heat  conduction 
equa  tion 


OT 

at 


=  a 


+  f  '4 

ox 


(4) 


The  justification  for  this  assumption  has  been  presented  in  several  papers  and 
need  not  be  present  'd  here.  However,  it  is  recognized  that  this  assumption  still 
may  not  be  valis.  More  discussion  on  this  point  will  be  presented  in  lattr  para¬ 
graphs  . 


At  the  solid  surface,  the  rate  at  which  energy  is  transferred  into  the 
solid  is  given  by 

-k(g)-F<Ts,  M)  (5) 

'x=+o 


TU  n  am  r  »-*  o  *- Vs  o  nf  f  v  f  v/vn  f*  V*  rw  one  »nKoqn  r*  (vnV »  «  O  — 

illC  lUuv.  uiwn  t,  v.Uii  tu  Atiw  t -»«w  v-  i  n.w  uo  w  *.  <.C  ^  *,  Un  fjww  f**“*>U  s.Oiiu/v*o 


tion  processes  well  as  the  energetics  of  any  surface  reaction/vaporization 
processes.  Also,  Ceep  within  the  p.opellant  the  temperature  is  not  influenced 
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by  the  combuscion  process.  Thus 

T  Tm  as  x  -►  ot  (6) 

Within  the  gas  phase,  ic  is  generally  possible  to  characterize  the  burning 
rate  by  the  general  expression 

mg  =  J(Ts,  P,  Tf)  •  o  (7) 

The  specific  details  of  the  gas-phase  combustion  will  determine  the  particular 
functional  form,  but  most  gas-phase  combustion  models,  be  they  diffusion  or 
premixed  flames,  will  follow  this  general  relationship. 

At  the  solid  surface,  the  rate  of  gasification  can  be  expressed  by 

r  =  H(T  ,  P)  =  m  A-  (8) 

s  s 

Equations  4  through  8  constitute  a  general  combustion  model  in  which  the 
only  assumptions  concerning  the  chemistry  and  physics  of  the  process  are  that 
the  bulk  solid  is  homogeneous  and  nonreactive.  These  equations  can  now  be  solved 
for  the  transient  burning  rate  and  gas  temperature  response  to  a  pressure  dis¬ 
turbance.  First,  since  m  =  m  ,  equations  7  and  8  can  be  combined  to  yield  an 

g  s 

expression  for  the  gas  temperature  in  terms  of  the  pressure  and  surface  tempera¬ 
ture 

Tf  =  G(Ts,  P)  (9) 

Then  the  gas  temperature  can  be  eliminated  from  equation  3  using  equation  9  to 
yield 


35 


UTC  2136-FR 


-kOT/^x)x=+o  =  F(Ts,  P)  (10) 

The  analysis  can  be  simplified  further  by  normalizing  the  characterizing  equa¬ 
tions  using  the  definitions 


e  =  (t-tJ/(t°  ~  t%)  r  =  rQ  x/a. 

2 

t  =  r  t/4 o  e  =  P/P 

o  o 

Hence,  equations  4,  8,  and  10  become 


+  H  (« 
n' 


u£ 

■y 


W0  =  */fco 


(11) 

(12) 

(13) 


-O0/c)Of=+o  =  F/Fo  -  W°  (14) 

and 

g  -*  o  as  t  -*  «> 

/ 1  o  \ 

It  should  be  noted  that  Hart  and  coworkers'  '  have  used  a  similar  approach 
for  characterizing  the  gas  phase  portion  of  the  combustion  zone.  However,  they 
used  an  Arrhenius  relation  to  characterize  the  surface  gasification  process. 

In  view  of  the  complexity  of  this  process,  especially  in  composite  propellants 
with  several  ingredients  decomposing  and  the  partially  molten  nature  of  the 
surface,  the  use  of  a  single  overall  pyrolysis  law  is  an  oversimplification  and 
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seems  unwarranted.  Also,  one  does  not  need  to  make  this  assumption  to  arrive 
at  a  useful  result,  as  the  following  development  will  show. 


To  solve  equations  12,  13,  and  14,  the  dependent  variables  can  be  expanded 
in  a  Taylor  series,  i.e., 


*  ’  1  +  *  '2  +  ^3  +  '  ' 


(15) 


In  a  similar  manner,  the  functions  H^,  and  can  be  expanded  by  a  Taylor 
series  to  yield,  for  example 

vll  2  o2H 

Hn(es,  )  =  l^e;  (-/)0+  '(^)o  +  £-  (t^)2(— y1)0 

1 0 

+  i2(“T1)o  +  2e2(5f)o  +  2  2  ^0+  r  ••• 


=  1+A.  +  Y~  •  P2+  • 


(16) 


Equations  15  and  16  can  then  be  substituted  into  equations  12,  13,  and  14 
to  yield  a  series  of  simultaneous  equations  characterizing  the  burning  rate 
behavior. 


2.2.2  Linear  Response 

The  linear  response  of  this  set  of  equations  can  be  determined  by  collect¬ 
ing  terms  of  the  order  \  i.e., 


*  =4*o 


(13) 
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It  should  be  uoted  that  a  detailed  knowledge  of  the  structure  of  the  gas 
phase  combustion  zone  or  the  surface  degradation  process  is  not  required  to 
characterize  the  acoustic  response  function.  The  combustion  parameters  repre¬ 
sent  the  partial  derivation  with  respect  to  surface  temperature  and  pressure 
of  both  the  burning  rate  and  the  energy  flux  into  the  solid.  These  general 
definitions  of  the  combustion  parameters  explain  why  Culick^^  was  able  to 
show  that  most  proposed  combustion  models  could  be  placed  in  the  form  of  equa¬ 
tion  25.  If  Q^=  0,  the  result  reduces  to  the  form  discussed  in  several  recent 
reports  and  papers,  with  a  slight  redefinition  of  the  parameter  B^. 

Further  manipulation  of  equations  19,  20,  and  21  shows  that  the  region  of 
intrinsic  stability  is  defined  by  the  relations 


4A!  >  Kj2 

(26) 

2AX  >  R1(R1-l) 

(27) 

If  these  conditions  are  not  satisfied,  any  perturbation  in  the  combustion  process, 
regardless  of  the  source,  will  cause  runaway  burning  rates.  They  form  a  boundary 
for  realistic  values  of  and  which  can  be  ascribed  to  real  propellants. 


2.2.3  Second  Order  Behavior 


Collecting  terms  of  the  order  A  yields 


(28) 
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+  0 


2 


(29) 


These  equations  are 
in  the  form  of  the  cross 
the  first  order  problem, 


linear  in  the  unknown  {?2>  nonlinearities  appear 
products  of  first  order  terms.  In  the  solution  of 
the  dimensionless  temperature  was  given  by 


0°  =  0^  exp(i.\r) 


However,  equations  28  and  29  were  derived  using  the  real  part  of  0“  .  Hence,  as 
Fried ly  has  shown,  the  cross  products  of  the  first  order  terms  must  be 
considered  as 


2  Re CV0.*)  (0  )2 

[Re  (0p]  =  - +  Re  ~ —  exp(2i\  t) 


(30) 


where  0 is  the  complement  of  0 


1' 


Hence  the  solution  of  the  dependent  parameters  in  equations  28  and  29  muse 
have  the  form 


sc 


Re  (6  )  =  Re(0j  +  Re  [9„  exp(2i>  :)] 

i.  c  t. 


(31) 
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Equations  30  and  31  show  that  the  solution  of  the  second  order  problem 
can  be  divided  into  two  parts  and  solved  separately.  One  part  reflects  the 
change  in  the  time  average  burning  rate  and  the  other  reflects  the  component 
which  oscillates  with  the  second  harmonic  frequency. 

After  making  the  appropriate  substitutions  into  equations  28  and  29,  the 
change  in  the  time  average  surface  temperature  and  hence  the  time  average 
burning  rate  can  be  determined. 

f/fQ  =  1  +  E(\)  .  x2  =  1  +  -  x2  Re (u^ /,  j2)  (32) 

The  detailed  expression  for  E(>)  is  presented  in  appendix  C.  It  should 
also  be  noted  that 

E  (0)  =  n(n-l)/4  (33) 


2.2.4  Third  Order  Behavior 

3 

Collecting  terms  of  the  order  A  yields 


j  b  6^ 


4  or 


v:2 


u0'j<9  dt? 

12,21.  c 

+  - u  — — 

v  o'  3  dr 


(34) 


(1?°)^  (C0)2.  1  r’2 D  r,3R 

/  a  , v  _  2ir  '•_r _ 1  _  .  ri  1  21  ,  ci  B3 

",c3  ^=+0  6  "  R3  +  2  D12  +  2  +  6 


+  <  ft  +  ei-J)DH+e;*2R2+  i'2B2+^3Rl  +  '3Bl 


(35) 
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The  proper  formulation  of  the  cross  products  in  equations  34  and  35  must 

(28) 

be  made  in  a  manner  similar  to  equation  30.  It  can  readily  be  shownv<"  '  that 

sc 

j  rn  ) 

Re[0^]  .Refe1^]  =  Re  j  +  — 2 —  1  exP  (^-^  (  +  terms  exP  (3i>.i)  (36) 

3  3(0 2  * 

Re(e®]  =  Re  [ — ^ -  (£!°)  exp(i>.i)j  +  terms  in  exp  (3i>.i)  (37) 

This  suggests  that  the  solution  for  should  take  the  form 

Re  [^3]  =  Re  [6”  exp(i\t)]  +  terms  in  exp(3i\t)  (38) 

Thus  the  solution  for  <r“,  and  the  third  order  burning  rate  response,  consists 
of  two  terms,  ore  which  oscillates  with  the  first  harmonic  and  one  which  oscil¬ 
lates  with  the  third.  The  first  term,  the  one  oscillating  at  the  first  har¬ 
monic,  is  of  particular  interest  because  it  represents  the  change  in  the  first 
order  response  with  pressure  amplitudes.  Thus 

Refo/cj)  =  ReC.^)  +  c*  F  (».)  (39) 

where  F(\)  represents  the  component  of  oscillating  with  exp(i\r).  The 
solution  of  equations  34  and  35  for  F(\)  is  presented  in  appendix  C. 

It  should  be  noted  that  the  principal  assumptions  in  the  preceeding  analysis 
are  that  the  solid  is  homogeneous  and  nonreacting  and  that  the  time  constants 
for  the  gas  phase  processes  are  sufficiently  fast  so  that  the  gas  phase  is  always 
in  equilibrium  with  the  instantaneous  pressure  and  solid  surface  conditions. 
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2.2.5  Utility  of  General  Mode 1 

Another  interesting  point  is  that  equations  20,  22,32,  and  39  can  be  used 
to  describe  the  linear  and  nonlinear  combustion  characteristics  of  any  combustion 
model  fitting  within  the  framework  of  equations  8,  12,  and  13.  The  simplified 
model  developed  early  in  this  program  may  be  used  as  an  example.  In  that 
analysis 

-k  —  =  K  +  2.  P  expC-E./R'I)  -  AH  .  f  (AO) 

x  4  4  v 

x=0 

(1)  (2)  (3) 

The  terms  on  the  right-hand  side  of  equation  40  represent 

(1)  =  heat  supplied  from  the  flame  zone  (gas -phase  reaction  zone 

a  finite  distance  away  from  the  surface) 

(2)  =  heat  generated  at  the  surface  by  heterogeneous  reactions 

between  the  solid  oxidizer  decomposition  gases  and  the 
exposed  fuel  matrix 

(3)  =  heat  released  because  of  gasification  of  the  propellant  at 

the  surface 

In  addition 

K  =  K  (P/P  )n 
o  o 

r  =  Z2  exp(-E2/RTs) 

For  this  particular  case 

A1  =  VTs'T«)/R(Ts)2 

Q,  =  0 

L 

R1  »  (  r :  - 

B  =  (n  ’’  +  r)/(1'«5) 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 
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Similarly,  one  could  wrice  the  equations  for  the  linear  and  nonlinear  behavior 
for  many  of  the  combustion  models  previously  reported  in  the  literature. 

2.3  COMPARISON  OF  THEORY  AND  EXPERIMENT 

A  comparison  of  the  experimental  results  described  in  section  2.1  was  made 
with  equations  20  and  22  for  two  reasons:  (1)  these  comparisons  yield  values 
of  the  linear  combustion  parameters  (A^,  B^,  and  R^)  ,  and  a  study  of  the 
effect  of  propellant  formulation  changes  on  these  parameters  will  provide 
additional  insight  into  the  combustion  process;  and  (2)  these  comparisons  pro¬ 
vide  some  idea  of  the  validity  of  the  basic  assumption  used  to  derive  equations 
20  and  22. 

2.3.1  Method  for  Deriving  Combustion  Parameters 

A  parameteric  study  of  equations  20  and  22  shows  th.Jt  the  maximum  value 
of  the  acoustic  response  function  depends  on  both  and  R^  and  that  the  fre¬ 
quency  at  which  this  maximum  occurs  depends  on  A^  and  Q^.  In  these  calculations, 
n  was  held  constant  and  was  determined  from  equation  22  for  the  particular 
set  of  Ap  R^ ,  Q^.  Since  them  is  no  simple  method  for  deriving  the  set  of 
combustion  parameters  graphically,  the  best  set  of  A^,  R^,  was  selected  by 
the  least-squares  approach  based  on  a  comparison  of  the  measured  and  predicted 
real  parts  of  the  response  function.  The  best  set  of  combustion  parameters 
derived  for  each  propellant  is  shown  in  tables  III  and  IV.  All  the  data  shown 
in  tables  III  and  IV  were  obtained  using  a  value  of  the  thermal  diffusivity  of 
3.0  x  10  4  in.^/sec.  The  effect  of  this  assumption  will  be  discussed  in  section 
2.3.5. 
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TABLE  III 


LINEAR  COMBUSTION  PARAMETERS 

AT  200 

PSIG 

Propellant  No. 

Formulation  No. 

A1  ±  2 

R1  ±  O' 

0 

Q1  +  0.2 

h 

1 

8501 

60 

7.0 

-0.3 

0.41 

2 

8502 

32 

4.4 

-0.3 

0.46 

3 

8525 

24 

5.2 

0.2 

0.48 

4 

8526 

26 

5.0 

-0.4 

0.45 

5 

8527 

No 

unique 

solution 

6 

8528 

14 

4.4 

1.0 

0.62 

7 

8529 

30 

4.6 

1.6 

18 

2.0 

3.6 

0.  57 

8 

8531 

12 

3.2 

2.4* 

0.95 

9 

8532 

36 

6.2 

1.3 

0.55 

10 

9164 

34 

3.6 

-0.5 

0.43 

11 

9165 

28 

5.2 

-0.1 

0.37 

12 

9167 

24 

1.8 

-0.3 

0.47 

13 

9168 

14 

3.0 

1.0 

0.57 

14 

9173 

18 

5.8 

0.6* 

0.47 

15 

9368 

38 

7.6 

0 . 6* 

0.35 

16 

9379 

36 

6.7 

-0.5* 

0.23 

17 

9395 

8 

3.4 

1 . 0* 

0.45 

18 

9715 

22 

5.4 

0.1 

0.31 

19 

9717 

10 

4.2 

0.1* 

0.29 

20 

9718 

18 

-0.4 

0.0 

0.53 

21 

9719 

12 

4.6 

0.3* 

0.45 

22 

9720 

28 

0.4 

-0.7* 

0.50 

23 

9754 

20 

6.2 

1.3* 

0.47 

*  Significant 

at  95%  confidence 

level . 

TABLE  IV 

EFFECT 

OF  PRESSURE  ON  COMBUSTION 

PARAMETERS 

(UTX-8501) 

Pressure  psia 

+  0.2 

Rx  +  0.2 

Qx  -  ±  0.2 

B, 

_ i 

100 

40 

7.8 

-0.7* 

0.27 

200 

60 

7.0 

-0.3 

0.41 

500 

20 

5.2 

0.0 

0.39 

*  Significant  at  95%  confidence  level. 
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Because  equations  20  and  22  are  nonlinear  in  the  combustion  parameters, 
the  process  of  selecting  the  minimum  value  of  the  error  estimate  (i.e.,  the 
sum  of  squares)  is  not  straightforward.  Of  the  several  methods  for  minimizing 
a  nonlinear  function,  the  grid  approach  was  selected  as  being  the  simplest  and 
most  straightforward.  The  region  of  ,  R^,  most  likely  to  contain  the  mini¬ 
mum  was  mapped,  and  the  set  of  A^,  R^,  yielding  the  minimum  error  estimate 
was  determined.  This  result  represents  a  "local”  minimum  which  may  not  be  the 
"universal"  minimum  which  represents  a  different  "best"  set.  This  is  illustrated 
by  the  fact  that  two  solutions  were  obtained  for  propellant  No.  7.  In  this  case, 
the  solution  A^  =  18,  R^  =  2.0,  =  3.6  yielded  the  lower  error  estimate. 

Propellant  No.  5  illustrates  another  problem  which  can  be  encountered. 

Only  four  measurements  were  obtained  for  this  propellant:  two  each  at  two 
different  frequencies.  Unfortunately,  there  was  not  a  significant  difference 
in  the  measured  response  functions  at  the  two  frequencies.  Thus,  one  could  draw 
the  response  function-frequency  curve  with  the  maximum  between  the  two  sets  of 
points,  or  one  could  consider  that  there  is  a  large  error  between  the  sets  and 
draw  the  curve  most  anywhere.  Therefore,  no  best  solution  could  be  found  for 
this  propellant. 

In  an  effort  to  ensure  that  the  "universal  best"  set  of  combustion  para¬ 
meters  was  determined,  a  wide  range  of  values  was  examined.  However,  there  is 
always  the  possibility  that  the  grid  was  too  large  and  that  points  of  low  error 
estimates  were  overlooked.  It  should  also  be  noted  that  propellants  where  the 
maximum  response  function  was  clearly  determined  did  not  present  any  problem, 
but  those  without  a  clearly  established  maximum  did  present  difficulties. 
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All  of  the  values  presented  in  tables  III  and  IV  satisfy  the  stability 
criteria  defined  by  equations  26  and  27  and,  thus , represent  values  which  are 
in  the  region  of  inherent  stability. 

2.3.2  Statistical  Significance  of 

Many  of  the  values  of  in  tables  III  and  IV  are  negative,  which  is  dif¬ 
ficult  to  explain  on  physical  ground.  However,  it  is  possible  that  the  values 
of  presented  in  tables  III  and  IV  are  not  significantly  different  than  zero. 

To  test  the  hypothesis,  the  data  were  examined  by  the  grid  approach  for 
the  case  of  =  0.  The  values  of  A^,  R^,  and  B.^  obtained  for  this  special 
case  are  shown  in  tables  V  and  VI.  The  error  estimates  for  the  general  case 
and  for  the  case  where  is  zero  were  compared  statistically  using  the  "F" 
test.  At  the  957*  confidence  level,  propellants  No.  8,  14,  15,  16,  17,  19,  21, 
22,  and  23(i.e.,  those  marked  with  a  flag  in  tables  III  and  IV)  have  values  of 
Q,  significantly  different  than  zero.  In  all  the  other  cases,  one  cannot  sepa¬ 
rate  the  derived  values  of  from  zero, 

2.3.3  Comparison  of  Data  Reduction  Method  with  Other  Approaches 

Several  alternate  data  reduction  approaches  to  the  one  used  in  this  program 
have  been  proposed  by  other  investigators.  With  the  exception  of  the  limited 

/ION 

efforts  of  Hart  and  coworkers,  '  these  comparisons  have  all  been  made  assuming 
-  0.  As  shown  in  section  2.3.2,  this  is  valid  in  some,  but  not  all,  cases. 

If  =  0,  then  one  approach  to  the  comparison  of  equation  20  with  experi¬ 
mental  data  involves  determining  both  the  real  and  the  imaginary  parts  of  the 
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TABLE  V 

LINEAR  COMBUSTION  PARAMETERS  AT  200  PSIG 
(Qi  =  0) 


TABLE  VI 

EFFECT  OF  PRESSURE  ON  COMBUSTION  PARAMETERS 
(UTX-8501  with  Qx  =  0) 
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» 

Propellant  No. 

Formulation  No. 

A1  ±  2 

RX  +  0.2 

h 

1  « 

1 

8501 

56 

6.6 

0.44 

2 

8502 

28 

4.2 

0.49 

• 

3 

8525 

24 

5.4 

0.43 

• 

4 

8526 

24 

4.6 

0.52 

5 

8527 

No 

unique  solution 

i 

6 

8528 

20 

5.2 

0.41 

7 

8529 

40 

6.8 

0.44 

i 

j  1 

8 

8531 

20 

5.6 

0.48 

. 

9 

8532 

52 

8.4 

0.36 

10 

9164 

30 

3.0 

0.48 

11 

9165 

26 

5.4 

0.37 

12 

9167 

20 

1.6 

0.49 

13 

9168 

28 

3.2 

0.46 

t 

14 

9173 

24 

6.0 

0.33 

15 

9368 

40 

8.2 

0.25 

2  * 

16 

9379 

32 

5.6 

0.32 

3  • 

17 

9395 

26 

6.2 

0.17 

1 

% 

18 

9715 

22 

5.6 

0.29 

] 

19 

9717 

12 

4.0 

0.29 

j  ’ 
i 

20 

9718 

18 

-0.4 

0.53 

1  * 

21 

9719 

14 

4.6 

0.38 

* 

22 

9720 

16 

1.6 

0.42 

23 

9754 

36 

7.8 

0.15 

Pressure,  psig 

A1  ±  2 

Rx  +  0.2 

h 

100 

44 

7.8 

0.42 

200 

56 

6.6 

0.44 

500 

20 

5.2 

0.39 

1 
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response  function  from  experimental  data.  For  the  T -burner, 
(°2)  (23) 

Fried lyv“  and  Becks tead  have  shown  in  earlier  work 


Strittma  ter, 


(20) 


it  ^-2— 
r  ~  4ocr 

o 


O  4rr , 

_£_d 


+  2jt  i 


(47) 


Since  the  growti:  and  decay  constants  (i.e. ,  r>  and  o^)  and  the  corresponding 
frequencies  are  determined  experimentally  for  each  point,  equating  the  real  and 
imaginary  parts  of  equations  25  and  47  provides  the  two  expressions  in  terms  of 
and 


A1  jRe^'^  ~  h^i*7  |  "  Rl  fRe(u/-')-n]  -  Im  (p.  /r ) 

!  e  j 

=  j  +  21, »(„/-:)  .J1  -  Re(u/«)  t  (>.xr)2  -  C'./)2] 


.  L 


-  Rx  [Re(ii/e) -n]  +  Im(u/<0 


jn,.^  -  2 Re (,*/>:)  \*  -  lm(u/0  -  (\11)2] 


(48) 


(49) 


Hence,  values  of  and  can  be  derived  for  each  data  point. 


This  approach  to  the  analysis  of  T-burner  data  has  a  serious  limitation 
resulting  from  experimental  uncertainties.  For  the  propellants,  shown  in  table 
IT, the  imaginary  part  of  equation  47  becomes 

Im(pA)  ~  400  (f  -f  )/f  (50) 

O 

Since  f  and  f^  have  similar  magnitudes,  the  experimental  errors  in  both 

and  are  magnified  in  the  calculation  of  Im(|i/r)»  This  results  in  significant 
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uncertainties  in  and  A1  when  they  are  derived  from  equations  48  and  49. 

Thus,  it  was  concluded  that  this  method  of  data  reduction,  though  more  rigorous, 
is  too  sensitive  to  small  errors  in  the  frequency  measurements.  For  this  reason, 
the  analysis  of  the  imaginary  part  of  the  response  function  was  not  included  in 
the  data  reduction  procedure. 

(23) 

In  an  effort  to  circumvent  this  difficulty,  Becks tead  and  Culick  have 
modified  this  approach.  Instead  of  using  the  imaginary  part  of  the  response 
function  in  equations  48  and  49,  they  assigned  Im(|i/e)  =  0.  To  determine  the 
effect  of  this  assumption,  some  of  the  data  obtained  in  this  program  were 
reduced  both  ways.  For  frequencies  below  800  to  1,000  cps,  the  two  methods 
yield  comparable  values  of  A^,  but  some  significant  changes  in  R^  were  noted. 

At  higher  frequencies  significant  differences  in  both  A^  and  R^  were  noted, 
thereby  casting  doubt  on  the  general  validity  of  the  modified  approach. 

(23) 

This  approximate  method  is  open  to  further  question  since  Becks tead 
found  that  the  linear  combustion  parameters  for  A-13  propellant  fell  outside 
the  region  of  intrinsic  stability.  This  result  is  clearly  not  satisfactory  on 
physical  grounds.  The  method  described  in  section  2.3.1  resulted  in  A^  =  16, 

R^  =  4.6,  and  =  0.5.  The  value  of  was  also  found  to  be  significantly 
different  than  zero.  These  values  are  well  within  the  region  of  intrinsic 
stability.  Even  if  the  results  were  based  on  =  0,  this  would  still  be  the 
case.  Therefore,  the  data  reduction  process  can  bias  the  derived  result  to 
the  extent  that  erroneous  conclusions  can  be  drawn. 
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(29) 

In  another  study,  Oberg  has  compared  experimental  data  with  equation 
20  using  only  the  real  part  of  the  response  function.  He  concluded  that  both 
the  T-burner  data  obtained  by  Horton^*^  and  the  L*  data  obtained  by  Beckstcad^^ 
can  be  adequately  correlated  by  equation  20. 


2.3.4  Approximate  Method  for  Data  Reduction 

In  the  course  of  studying  various  data  reduction  methods,  a  quick  and 
accurate  method  was  developed  in  which  =  0  is  assumed.  Under  these  conditions, 
parameteric  solutions  of  equation  25  show  that,  to  a  good  approximation,  the 
response  is  a  maximum  when  Im(p^/e^)  =  0.  Hence, 

A1  =  l!  4(>*ir  '  l)  (51) 

where  >.^  is  evaluated  at  the  frequency  of  maximum  response. 

A  plot  of  equation  51  is  shown  in  figure  20  where  R,C  is  the  critical  value  of 
R1  for  intrinsic  stability  determined  ii.  equation  27,  Therefore,  the  parameter 
A^  can  be  determined  directly  from  experimental  data  by  noting  the  dimensionless 
frequency  at  which  Re(^/c)  becomes  a  maximum.  The  parameter  R^  can  then  be 
determined  from  the  maximum  response  function  and  burning  rate  pressure  exponent 
by  substituting  equation  51  into  equations  20  and  22,  yielding 


1 

n 


Re(u/<r) 


max 


(l-R^+Ap/k^ 


-R^  +  4Aj\^ 


(52) 


Equation  52  is  shown  graphically  in  figure  21. 


The  values  derived  by  this  method  are  nearly  identical  to  those  presented 
in  tables  V  and  VI  which  were  derived  by  the  least-squares  approach. 
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2.3.5  Analysis  of  Error  Effects 

For  reasons  of  simplicity  and  economy,  the  approximate  method  presented 
in  section  2. 3. A  vas  used  in  this  part  of  the  study.  The  effect  of  errors  in 
the  dimensionless  frequency  at  the  maximum  response  (i.e.,  burning  rate,  dif- 
fusivity,  and  dimensional  frequency  errors)  were  determined  by  an  arbitrary  16% 
reduction  in  a  and,  therefore,  the  dimensionless  frequency  at  the  maximum  re¬ 
sponse.  The  results  show  that  an  error  of  16%  in  X  results  in  a  16%  error  in 
and  an  8%  error  in  R^. 

Arbitrary  changes  in  ReCp/e)  and  n  were  evaluated  and  the  resulting  effect 
on  R^  is  shown  in  table  VII.  Note  that  A^  is  determined  by  X  alone  and  is  not 
subject  to  these  errors. 


TABLE  VII 

EFFECT  OF  n  AND  Re(p/e)  ON  R, 

max  i 


5.5 

5.7 

5.9 

0.49 

|  6.60 

6.77 

6.93 

0.55 

Rx  1  6.00 

6.19 

6.37 

0.60 

I  5.43 

5.70 

5.90 

The  results  show  that  errors  in  the  experimental  measurements,  i.e,,  X,  n,  and 
Re(p/e)  ,  result  in  approximately  the  same  percentage  error  in  the  derived 
combustion  parameters.  Hence,  the  procedure  for  deriving  the  combustion  para¬ 
meters  does  not  significantly  amplify  or  attenuate  experimental  uncertainties. 
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The  response  function- frequency  curve  was  calculated  for  propellant  No.  1 

using  the  derived  combustion  parameters  shown  in  table  IV.  Figure  22  shows 

that  the  shape  of  the  predicted  response- frequency  curve  is  not  significantly 

altered  by  an  error  in  a.  Figure  23  shows  the  effect  of  n  and  Re(nA)  on 

m8X 

Ih..:  predicted  response  function  frequency  curve.  The  pressure  exponent,  n, 
influences  the  shape  of  the  curve  at  the  frequency  extremes  but  has  little  in¬ 
fluence  near  the  maximum;  Re(u/Omax  ^as  roore  influence  on  the  peak  region  of 
the  curve  but  little  influence  on  the  frequency  extremes. 

Figures  22  and  23  show  that  the  response  function  derived  from  equations 
20  and  22  adequately  predicts  the  observed  shape  of  the  frequency-response 
function  curve.  Other  comparisons  between  the  data  on  the  curves  predicted  by 
equations  20  and  22  using  the  combustion  parameters  shown  in  table  III  are  in 
agreement.  This  leads  to  the  conclusion  that  equations  20  and  22  adequately 
represent  the  linear  combustion  behavior  of  the  propellants  listed  in  table  II. 

It  would  be  interesting  to  compare  equations  20  and  22  with  the  extensive 
data  obtained  at  NWC,  but  these  data  have  not  been  publicly  distributed. 

2.3.6  Interpretation  of  Combustion  Parameters 

It  is  interesting  to  speculate  on  the  physical  significance  of  the  various 
linear  combustion  parameters  and  the  numerical  values  of  the  associated  physical 
parameters.  The  parameters  associated  with  the  burning  rate  of  the  solid  phase, 
i.e.,  the  function  H(Ts,P),  can  be  considered.  If  one  supposes  that 

is  =  H(Ts,P)  =  Z2  P^l  exp(-E2/RTs)  (53) 
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(a  form  which  has  been  proposed  in  almost  all  combustion  models  but  which  has 
not  been  convincingly  justified),  then 

Ax  =  E2(T“-To)/RT“s2  (54) 

BL  =  n1  (55) 

Given  equation  53,  the  parameter  A^  is  a  measure  of  the  surface  temperature 
and  actuation  energy  of  the  reaction.  Assuming  an  initial  propellant  tempera¬ 
ture  of  300°K  and  a  steady-state  surface  temperature  of  750°K,  a  value  of  A^ 
equal  to  30  corresponds  to  an  E£  of  74.5  Kca 1/mole.  If  the  surface  temperature 
is  assumed  to  be  850°K,  than  an  A^  of  30  corresponds  to  an  of  78.1  Kcal/mole. 
Thus,  in  the  case  of  equation  53,  the  variations  in  A^  shown  in  table  III  and  IV 
reflect  formulation  effects  on  E,  versus  the  effects  on  the  surface  temperature. 
At  a  750°K  surface  temperature,  the  activation  energies  derived  from  table  I 
range  from  20  Kcal/mole  for  propellant  No.  17  to  145  Kcal/mole  for  propellant 
No.  1.  These  are  extremely  high,  especially  in  view  of  the  values  reported  for 

the  decomposition  and  gasification  of  various  ingredients.  The  results,  though 

(23) 

consistent  with  nonacoustic  combustion  instability  data,  suggest  that  if 
equation  53  applies  there  is  a  complex  series  of  surface  reactions  and  that  the 
formulation  changes  shift  the  balance  rather  drastically. 

Further  evidence  of  this  shift  in  the  reaction  process  is  provided  by  the 
catalyst  data  in  table  III.  These  data  show  the  addition  of  catalysts  to  the 
formulation  decreases  the  value  of  A^;  the  higher  the  catalyst  concentration 
the  lower  the  value  of  A^.  At  the  same  time,  is  increased  from  -0.5  to 
approximately  0.1  when  the  catalysts  are  present  in  concentrations  greater  than 
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0.25%.  Since  Che  parameters  and  reflect  the  solid  phase  contribution  to 
the  mass  burning  rate,  the  effect  of  burning  rata  catalysts  influence  this  part 
of  the  combustion  process  in  addition  to  the  gas  phase  reaction  effects . 

The  data  in  table  III  indicate  that  the.  effect  of  pressure  is  significant 
at  the  higher  solids  loading.  At  78%  oxidizer  loading,  only  two  of  the  propel¬ 
lants  showed  values  significantly  different  than  zero.  At  84%,  however, 
only  two  of  the  nine  propellants  had  values  of  zero.  Thus,  the  coatings 
had  little  effect  on  the  pressure  dependence  of  the  solid  phase  aspects  of  the 
mass  burning  rate. 

A  perplexing  aspect  of  these  data  ir  the  fact  that  some  of  the  significant 
nonzero  values  of  are  negative,  while  they  are  greater  than  one  in  other  ■ 
cases.  The  latter  result  may  be  somewhat  surprising,  but  it  can  be  explained 
in  terms  of  possible  chemical  reaction  mechanisms.  The  explanation  of  the 

former  result  (i.e.,  that  Q^'O)  is  not  possible  in  this  manner.  This  difficulty 

in  explaining  the  activation  energy  and  the  order  of  the  reaction  leads  one  to 
suspect  the  validity  of  equation  53  in  characterizing  propellant  combustion. 
However,  logical  alternates  to  equation  53  are  equally  difficult  to  propose. 

As  further  evidence  of  the  complexity  of  physically  interpreting  equation 
8,  consider  the  results  presented  in  cable  TV,  The  effect  of  pressure  on 
suggests  that  as  the  pressure  is  increased,  0^  becomes  less  negative  (therefore 
more  positive).  At  first,  one  would  ascribe  this  change  to  the  increasing  rate 
of  the  self  combustion  of  the  ammonium  perchlorate  at  the  higher  pressure.  The 

effect  of  pressure  on  A^  also  supports  the  concept  of  a  shifting  reaction 
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mechanism  with  pressure.  However,  one  is  still  faced  with  the  difficulty  of 
explaining  negative  values  of  at  the  lower  pressures  and  the  parabolic 
effect  of  pressure  on  A^. 


The  situation  regarding  and  is  even  more  complex.  Tc  demonstrate 
this  point,  the  first  order  Taylor  expression  of  equation  5  is 


Equation  56  can  then  be  combined  with  the.  first  order  expansions  of  equations 
7  and  8  to  provide  alternate  definitions  of  R^  and  in  terms  of  the  functions 
K(Ts,  P,  Tf),  H(Ts,  P),  and  J(Tg,  P,  Tf) . 


(oK/oTs  )  q  (<5  J  /bT  f  )  o  +  (bK/^Tf)Q[0H/?>Ts)o-(aj/aT8)o]; 

i  o  C  (bJ/bT  ) 
o  p  r  o 


B  =  P 
1  o 


j(bK/bP)o(bJ/bTf)Q  +  (bK/bTf)ot(bH/dP)o-(bJ/bP)o] 


i  o  C  (T  -T  )(dJ/dT.) 
o  p  s  00  to 


Considering  the  wide  range  of  proposed  definitions  for  the  parameters  in 
equation  56  one  finds  that  generally  (bH/bT  )>(bJ/bT  ),  that  (bK/bT-)>0,  and 
that  J/bT^)>0.  Hence,  one  can  conclude  that  if  R^'O,  then  (bK/bTs)<0.  This 
would  be  the  case,  for  example,  if  equation  53  applies  and  the  process  is  endo¬ 


thermic.  However,  since  the  experimental  value  of  Rp>0,  this  implies  that 
approximately  (oK/bT  )>0.  This  then  suggests  that  the  surface  processes  are 

S  —  /  v  j  \ 


at  least  nearly  thermally  neutral 


,.  bK  (m\  CM  bJ  w/bJ\. 

U.e.,  £Xg  ~  l  oTfyl  0Tg  '  bTg  J/  lbTf  J 
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or  perhaps  even  somewhat  exothermic.  The  magnitude  of  the  various  partial 
derivatives  determine  the  extent  to  which  this  is  true.  However,  it  would 
appear  that  the  values  of  are  reasonable  in  terms  of  the  various  proposed 
definitions  from  the  several  available  combustion  models.  Similarly,  the 
values  of  derived  from  the  data  are  also  quite  reasonable.  Either  the  granu¬ 
lar  diffusion  visible  or  the  laminar  premixed  flame  model  predicts  pressure 
indices  between  0  and  1,  which  is  consistent  with  the  observations. 

Examination  of  the  effects  of  the  various  formulation  changes  on  R^  and 
B^  yields  some  interesting  results.  The  effect  of  coatings  on  the  solid  oxidi¬ 
zer  can  be  considered  first.  The  coatings  reduced  the  temperature  sensitivity 
of  the  heat  transfer  into  the  solid  phase.  Viton  A  was  the  least  effective,  and 
Kel-F,  Hypalon,  and  ethyl  cellulose  all  had  nearly  the  same  effect.  It  is,  of 
course,  impossible  to  determine  whether  the  effect  is  the  result  of  a  reduction 
in  possible  exothermic  surface  reactions  (the  original  intent  of  the  experiment) 
or  the  result  of  decreased  mixing  within  the  gas  phase  reaction  zone.  There 
appears  to  be  little  effect  of  the  coatings  on  B^,  so  no  clues  can  be  obtained 
from  studying  those  results.  It  should  also  be  noted  that  the  coatings  were 
porous,  thereby  lessening  the  resulting  effects  in  comparison  to  a  nonporous 
coating. 

The  influence  of  the  catalysts  is  also  interesting  to  consider.  In  all 
cases,  the  values  of  R^  were  substantially  reduced  by  the  addition  of  the 
catalysts.  When  combined  with  their  influence  or  the  activation  energy  of  the 
surface  gasification  process,  the  data  suggest  the  effect  is  to  reduce  the 
activation  energies  of  the  various  reactions.  An  additional  influence,  which 
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effects  R^  but  not  Ap  results  from  increasing  the  gas  phase  damping.  Both 
catalysts  exist  as  solids  within  the  gas  phase  combustion  zone  and  could  behave 
in  much  the  same  manner  as  aluminum. 

Another  aspect  of  the  catalyst  data  is  that  the  addition  of  copper  chromite 
results  yield  substantially  lower  values  of  R-  and  higher  values  of  B. .  It  may 

JL  JL 

be  that  copper  chromite  influences  different  processes  characterized  in  and 
B^  than  iron  oxide. 

The  influence  of  pressure  on  R^  and  B^,  at  least  for  propellant  No.  1, 
seems  to  be  counter  to  expectations.  Increasing  the  pressure  appears  to  de¬ 
crease  R^  in  the  region  where  one  would  expect  the  self-deflagra tion  of  AP  to 
increase  R^.  However,  the  shift  in  the  complex  reaction  path  and  the  resulting 
change  in  activation  energies  could  cause  the  effect  to  be  opposite  to  the 
anticipated  effect. 

2.3.7  Nonlinear  Behavior 

The  nonlinear  behavior  of  the  combustion  model  has  also  been  examined.  The 
original  approach  to  the  study  of  the  nonlinear  constants  was  to  attempt  their 
derivation  by  extrapolating  the  linear  results.  However,  the  lack  of  decisive 
linear  results  made  '*  pproach  impractical.  Therefore,  two  alternatives  were 
available:  (1)  study.  ig  the  effects  of  the  nonlinear  constants  separately;  (2) 
studying  the  nonlinearities  within  the  context  of  one  particular  combustion  model. 
Both  approaches  have  been  attempted,  and  the  results  from  each  approach  are 
presented  in  the  following  paragraphs. 
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2.3. 7.1  Parameteric  Studies  of  General  Analysis 

Equations  20,  22,32,  and  39  show  that  a  large  number  of  constants  are 
required  to  characterize  the  complete  behavior.  However,  a  parameteric  study 
of  the  shift  in  time  average  burning  rate  is  straightforward  if  one  keeps  the 
linear  behavior  fixed.  Equation  32  shows  that  if  -  0,  which  is  reasonable 
in  the  T-burner,  then  only  two  second  order  constants,  X^  and  W^,  are  required. 
Calculations  were  then  made  assuming  =  20,  n  =  0.40,  =  4.0,  and  =  0. 

The  effect  of  and  is  shown  in  figure  24. 

These  results  show  that  both  increases  and  decreases  in  the  time  average 
burning  rate  can  be  explained  within  this  analytical  framework  depending  on  the 
nonlinear  constants.  It  should  be  noted  that  in  both  cases  the  zero  frequency 
corresponds  to  the  shift  consistent  with  the  equation  33.  Only  limited  experi- 

(311 

mental  data  are  available  with  which  to  compare  these  predicted  results.  Eisel 
has  found  that  -0.5  E(\)^  0  for  both  composite  and  double-base  propellants. 

When  compared  with  figure  24  and  other  parameteric  results  not  shown,  the  data 
indicate  that  is  less  than  200  for  these  propellants. 

Turning  to  third  order  effects,  parameteric  studies  show  that  F(\)  can 
have  a  wide  range  of  values  depending  on  both  second  and  third  order  constants. 

The  effects  are  too  numerous  for  a  meaningful  study  of  the  behavior. 

2. 3. 7. 2  Parametric  Study  for  a  Particular  Combustion  Model 

The  model  selected  to  study  the  nonlinear  combustion  behavior  was  the 
simplified  model  presented  in  section  2.2.2.  This  model  does  have  its  limita¬ 
tions,  particularly  in  the  analysis  of  the  transient  flame  temperature.  However, 
the  derivation  of  equation  25  shows  that  this  effect  can  be  treated  separately 


63 


UTC  2136-FR 

from  the  analysis  of  the  transient  burning  rate.  Therefore,  the  results  of 
this  nonlinear  analysis  are  somewhat  more  general  than  the  combustion  model 
sugges  ts . 

An  examination  of  the  higher  order  equations  28  through  39  reveals  that 
six  combustion  parameters  are  required  to  characterize  the  nonlinear  behavior 
of  the  simplified  model  and,  as  long  as  they  represent  an  independent  set,  any 
combination  of  six  is  sufficient.  In  the  present  studies,  the  three  parameters 
characterizing  the  linear  behavior  were  retained.  The  three  remaining  para¬ 
meters  are  the  surface  temperature  ratio  ({J  =  T^  /T“  ),  the  pressure  exponent 
of  the  gas  phase  heat  flux  (n1),  and  the  ratio  of  activation  energy  of  exothermic 
surface  reaction  to  that  of  the  overall  pyrolysis  reaction  at  the  surface 
(  A  =  E4/E2).  These  three  additional  parameters  represent  physical  characteris¬ 
tics  of  the  combustion  process  which  can  be  estimated  within  reasonable  bounds. 

For  example,  it  is  generally  believed  that  typical  surface  temperatures  are  of 

(32) 

the  order  of  750°K,  '  thereby  providing  a  reasonable  estimate  for  0.  Likewise, 

n*  must  lie  between  zero  and  one.  A  is  the  least  \nown  parameter,  but  in  the 
light  of  the  high  values  of  A  should  also  lie  between  zero  and  one, 

A  parameteric  study  of  the  higher  order  effect  has  been  made  using  this  set 
of  combustion  parameters.  The  results  are  shown  in  figures  25  through  27  for  the 
case  where  the  linear  response  (acoustic  response  function)  corresponds  to  propel¬ 
lant  No.  1  in  table  II.  Figures  25  and  26  show  E(\)  is  always  negative  and  follows 
the  inverse  frequency-acoustic  response  function  pattern  (i.e,,  E(\)  is  small 
at  lov,  and  high  frequencies)  and  reaches  a  minimum  at  the  same  frequency  in 
which  the  response  function  is  a  maximum.  The  results  show  the  minimum  value 
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of  E(>.)  becomes  increasingly  negative  for  decreasing  values  of  the  gas  phase 
pressure  exponent  (n1),  as  shown  in  figure  25.  The  effect  of  A  on  E(X.), 
which  is  shown  in  figure  26, is  more  complex.  For  a  >0.25  the  minimum  value  of 
E(\)  increases  with  increasing  values  of  the  activation  energy  ratio  A  .  For 
A  -'0.25,  the  minimum  value  of  E(\)  increases  with  decreasing  values  of  A  . 
Other  calculations  show  that  the  temperature  ratio,  0,  has  little  effect  on 
E(\). 


These  results  can  be  explained  in  terms  of  the  physics  of  the  combustion 
model.  If  n'  is  less  than  unity,  the  energy  transfer  from  the  gas  phase  will 
not  be  symmetric  about  the  mean  value  during  one  pressure  cycle.  In  fact, 
there  will  be  a  net  decrease  in  energy  transfer,  depending  on  n',  which  results 
in  the  decrease  in  steady  burning  rate.  The  effect  of  A  results  from  the 
asymmetry  of  the  surface  heat  liberation  term  as  a  consequence  of  the  surface 
temperature  contribution.  If  A  is  significantly  less  than  unity,  the  increase 
in  surface  temperature  does  not  increase  the  burning  rate  and  energy  liberation 
in  direct  proportion.  The  effect  of  A  is  complicated  by  the  fact  that  the  heat 
release  from  surface  reactions  depends  on  both  surface  temperature  and  pressure. 
This  dual  dependence  accounts  for  the  nonmonotonic  effect  of  A  .  For  A  =  0, 
the  surface  reaction  effect  depends  only  on  pressure  to  the  first  power  and  it 
oscillates  symmetrically  about  the  mean.  This  symmetric  effect  counters  the 
asymmetric  behavior  of  the  gas  phase  flux  and  reduces  the  overall  asymmetry  of 
the  energy  transfer  to  the  solid. 

The  behavior  at  the  frequency  extremes  can  also  be  explained  on  this  basis. 
At  low  frequencies,  asymmetry  results  because  the  burning  rate  pressure  exponent 
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is  less  than  one.  At  the  high  frequency,  asymmetry  results  from  the  behavior 
of  the  gas  phase  flux  alone,  since  the  surface  temperature  can  not  follow  the 
pressure  oscillations  (i.e.,  Re(p/?)-»0  as  V>  oo) .  This  contrasts  the  high- 

/og\ 

frequency  behavior  reported  by  Friedly.  ' 

Similar  calculations  were  made  for  the  corresponding  PBAN/AP  propellant 
(i.e.,  propellant  No.  13).  The  results  show  similar  behavior  to  that  shown 
in  figures  25  and  26,  except  that  the  magnitude  of  the  effects  are  approximately 
half  that  shown  in  these  figures.  Therefore,  the  parameters  which  characterize 
the  linear  behavior,  particularly  R^,  have  a  significant  influence  on  the  non¬ 
linear  behavior  as  well. 

The  maximum  change  in  E(\),  shown  in  figures  25  and  26,  represents  a  small 

shift  in  burning  rate.  Typical  values  for  in  the  T-burner  are  0.1  to  0.3,  and 

the  maximum  decrease  in  burning  rate  is  approximately  2%.  This  value  is  of  the 

(31) 

same  order  of  magnitude  as  the  experimental  observations  of  Eisel.  Further¬ 

more,  the  change  in  burning  rate  is  shown  to  be  quadratic  in  pressure  amplitude, 
which  also  agrees  with  experimental  results.  Unfortunately,  the  accuracy  of  the 
data  is  not  sufficient  to  make  meaningful  quantitative  comparisons  with  theory. 

Another  aspect  of  the  nonlinear  studies  concerns  the  changes  which  occur 
in  the  acoustic  response  function  with  increasing  pressure  amplitudes.  The 
effects  of  0,  n‘,  and  A  on  F(\)  were  determined  through  equation  23,  using  the 
values  of  R^,  A^,  and  n  for  propellant  No.  1  (thereby  matching  the  experimentally 
determined  linear  response).  Figure  27  shows  th',  effect  of  A  on  F(\).  The 
effects  of  n*  and  $  are  found  to  be  essentially  negligible.  It  should  be  noted 
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thac  the  curves  for  0  <  A  '0,50  cannot  be  distinguished  from  the  curve  for 
A  =  0.50.  The  results  also  show  that  F(\)  approaches  zero  at  both  high  and 
low  frequency.  The  high-frequency  behavior  results  from  the  inability  of  the 
solid  to  respond  to  the  pressure  oscillations  while  the  low-frequency  limit 
indicates  that  the  burning  rate  pressure  exponent  n  is  independent  of  pressure. 

(33) 

The  results  shown  in  figure  27  are  similar  to  Friedly's'  '  predictions, 
although  there  are  some  important  differences.  F(?v)  is  negative  at  frequencies 
which  are  less  than  the  frequency  at  which  the  maximum  response  occurs.  Friedly's 
results  would  lead  to  the  opposite  conclusion.  One  possible  explanation  is  that 
different  methods  were  used  to  solve  the  differential  equations  and  calculate 
F(\):  Friedly  used  the  approximate  method  of  discretization,  in  contrast  to  the 
exact  variation  of  parameters  approach  used  in  this  study. 

The  large  negative  values  in  the  intermediate  pressure  region  show  that  the 
combustion  process  is  self-limiting  in  the  T-burner  configuration.  As  the  pres¬ 
sure  amplitude  increases  because  of  the  linear  behavior  in  the  T-burner,  the 
response  function  decreases.  This  reduces  the  energy  input  driving  the  acoustic 
pressure  and  eventually  results  in  an  upper  limit  in  the  acoustic  pressure.  By 
setting  Re(^/c)  equal  to  zero  m  equation  39,  these  limiting  pressure  amplitudes 
have  been  calculated  for  the  c  see  shown  in  figure  27.  The  results,  which  are 
shown  in  figure  28  reveal  that  limiting  amplitudes  of  0.4  would  be  expected  on 
the  basis  of  combustion  oehavior  alone  at  typical  T-burner  frequencies.  Since 
these  calculations  do  not  consider  the  effect  of  acoustic  losses,  the  predictions 
are  reasonable  in  the  light  of  experimentally  observed  limit  amplitudes.  Further¬ 
more,  calculations  show  the  effect  of  n'  and  0  to  be  insignificant,  and  the  effect 
of  A  cn  the  limit  amplitude  can  be  anitiepated  from  figure  27. 
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2.3.8  Effect  of  Pressure 

The  effect  of  pressure  on  the  transient  behavior  is  also  of  interest  in 

(34) 

light  of  the  data  of  Ibiricu.  This  effect  can  only  be  studied  using  a 

specific  combustion  model.  This  was  done  with  the  model  described  by  equations 
40,  41,  and  42  because  of  the  simplicity  and  economy  of  calculations.  This 
combustion  model  predicts  that  the  response  function  increases  with  pressure 
and  that  the  frequency  at  which  the  maximum  response  occurs  also  increases 
without  limit.  Using  the  previously  reported  method  to  correct  for  changes 
in  propellant  surface  temperature  with  pressure*^  the  predicted  result  of 
increasing  pressure  is  shown  in  figure  29  for  propellant  No.  1.  In  this  cal¬ 
culation  A  =  0.25,  n*  =  0.40,  and  0  =  0.40,  while  the  values  of  R^,  A^,  and 

n  at  200  psi  are  those  shown  in  table  II.  This  effect  is  consistent  with  the 

(34) 

stability  behavior  observed  by  Ibiricu  at  the  lower  pressures.  However, 
it  the  higher  pressure,  Ibiricu  observed  the  response  function  to  pass  through 
a  maximum  and  then  to  decrease  with  further  increases  in  pressure.  The  linear 
behavior  shown  in  figure  29  shows  no  such  maximum. 

Further  study  of  the  pressure  influence  on  the  nonlinear  behavior  reveals 
a  partial  explanation  of  this  dichotomy.  The  effect  of  pressure  on  the  limit¬ 
ing  pressure  amplitudes  for  propellant  No.  1  using  the  values  of  A  ,  n',  and 
0  is  shown  in  figure  30.  From  this  plot,  it  is  apparent  that  increasing  the 
pressure  sharply  reduces  the  limiting  pressure  amplitude. 
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These  results  now  provide  a  possible  explanation  for  the  data  of  Ibiricu. 

At  the  lower  pressure  (i.e.,  approximately  200  psia)  the  limiting  pressure 
amplitudes  are  sufficiently  large  to  permit  meaningful  measurement  of  the  ampli¬ 
tude  growth  constants  in  the  linear  region.  Therefore,  the  reported  response 
functions  are  truly  characteristics  of  the  propellant  behavior,  In  this  area 
Ibiricu's  data  and  the  data  obtained  for  propellant  No.  1  in  the  course  of  this 
study  are  consistent  with  the  predicted  linear  behavior. 

As  the  pressure  is  increased,  however,  the  limiting  acoustic  pressure 
amplitude  decreases.  Considering  the  fact  that  nonlinear  effects  become  signifi¬ 
cant  when  the  pressure  amplitude  reaches  50%  of  the  limiting  value,  the  region 
in  which  meaningful  response  function  measurements  can  be  made  is  sharply  reduced 
with  increasing  pressure.  Furthermore,  if  nonlinear  effects  are  mistakenly  in¬ 
cluded  in  the  linear  interpretation  of  the  data,  equation  39  shows  that  the 

resulting  response  function  will  be  less  than  the  true  value.  Hence,  it  is 

(34) 

possible  that  the  measurements  reported  by  Ibiricu  might  reflect  the  apparent 
response  function  and  not  the  true  response  function. 


To  determine  if  the  magnitude  of  this  effect  is  sufficient  to  account  for 
the  data,  equations  2  and  27  have  combined  to  yield  the  following  expression  for 
the  growth  of  the  acoustic  pressure  amplitude  in  the  T-burner 


Experimental  results  from  T-burner  measurements  have  indicated  that,  as  a  first 

3/4  1/4  (35) 

approximation,  o  =  Kf^  /(PQ).  '  An  integration  of  equation  59  yields 
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the  amplitude  of  the  acoustic  pressure  as  a  function  of  time.  Values  of  the 
apparent  response  function  were  derived  from  these  predicted  curves  by  measuring 
the  apparent  growth  constant  at  a  pressure  amplitude  of  5  psi.  The  difference 
between  the  resulting  apparent  response  function  and  the  true  value  (determined 
at  zero  pressure  amplitude)  was  negligible  over  a  wide  pressure  region. 

This  calculation  is  somewhat  inconclusive,  since  the  effect  of  pressure 
on  the  combustion  parameters  was  determined  by  the  limited  specific  combustion 
model.  Thus,  one  can  conclude  that  this  specific  model  does  provide  the  magni¬ 
tude  of  the  pressure  effect  required  to  account  for  Ibiricu's  observations.  How¬ 
ever,  other  combustion  models  which  consider  the  transient  flame  temperature 
behavior  might  account  for  the  experimental  observations. 
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3.0  CONCLUSIONS 

A.  Transient  combustion  of  composite  solid  propellants  can  be  character¬ 
ized  by  assuming  the  transient  behavior  of  the  solid  phase  to  be  that 
of  a  homogeneous  nonreacting  solid. 

B.  The  gasification  of  the  solid  phase  was  found  to  be  an  extremely  complex 
process.  For  all  propellants  tested,  this  process  was  found  to  be 
strongly  dependent  upon  surface  temperatures,  more  so  than  could  be 
accounted  for  by  an  Arrhenius  expression  with  a  reasonable  activation 
energy.  For  propellants  having  an  847.  loading  of  AP,  the  gasification 
process  was  also  dependent  on  pressure  and  surface  temperature.  Some 
propellants  showed  the  gasification  increasing  with  pressure,  while 
others  showed  a  decrease. 

C.  The  rate  of  energy  transfer  into  the  solid  was  found  to  increase  with 
increasing  surface  temperature  and  increasing  pressure.  The  tempera¬ 
ture  effect  is  characteristic  of  a  thermally  neutral  or  slightly 
exothermic  surface.  The  pressure  effect  is  consistent  with  that  pre¬ 
dicted  by  laminar  premixed  and  granular  diffusion  flame  models. 

Coating  the  solid  oxidizer  reduced  the  temperature  sensitivity  of  the 
energy  transfer,  with  the  magnitude  of  the  effect  depending  on  the 
coating.  However,  coatings  had  little  effect  on  the  pressure  sensi¬ 
tivity. 


78 


UTC  2136-FR 


D.  Parameteric  studies  were  conducted  in  which  the  linear  behavior  was 
held  constant  and  the  effects  of  the  second  and  third  order  coef¬ 
ficients  were  determined.  The  results  show  that  both  positive  and 
negative  values  of  E(A)  can  be  obtained,  depending  on  the  magnitude 
of  the  constants.  Since  Eisel  has  observed  E(A)  to  be  negative  in 
his  studies ,  this  defines  the  region  of  reasonable  second  order  coef¬ 
ficients.  Additional  studies  show  that  F(A)  is  negative  for  reason¬ 
able  second  and  third  order  coefficients.  Therefore,  the  combustion 
process  itself  can  limit  the  acoustic  pressure  amplitude  in  the  T- 
burner.  On  the  basis  of  these  calculations,  limiting  amplitudes  of 
0.1  to  0.5  can  be  expected  without  consideration  of  acoustic  losses. 
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APPENDIX  A.  ADDITIONAL  EXPERIMENTAL  STUDIES 

1.0  PHOTOGRAPHIC  STUDIES 

A  sample  holder  was  constructed  with  viewing  windows,  as  shown  schemati¬ 
cally  in  figure  A-l.  The  1. 5-in. -diameter  propellant  disk  was  located  in 
the  center,  and  a  purge  sys  ten.  was  located  between  the  sample  and  the  windows. 
Tests  were  conducted  using  both  Vycor  and  Plexiglas  windows  with  and  without 
the  purge  system.  The  windows  were  sealed  to  their  mounting  with  gaskets.  A 
photograph  of  the  final  assembly  is  shown  in  figure  A-2. 

The  camera  used  in  these  studies  was  a  Hycan  16-mm  rotary  prism  camera 
with  a  framing  rate  capacity  to  4,000  frames/sec,  which  will  provide  several 
frames  during  a  complete  cycle  at  the  lower  frequencies.  The  initial  photo¬ 
graphs  were  taken  at  a  framing  rate  of  500  frames /sec.  A  quartz -iodine  lamp 
was  used  for  backlighting  the  viewing  section. 

The  greatest  problem  encountered  in  the  photographic  studies  was  the 
rapid  buildup  of  carbonaceous  material  on  the  windows.  To  eliminate  this  prob¬ 
lem,  several  changes  were  made  in  the  experiments.  A  clean-burning  propellant 
consisting  of  84'  AP  was  selected.  The  burning  rate  at  200  psig  was  0.130 
in. /sec.  The  windows  were  also  recessed  from  the  tube  walls  to  allow  the  nitro¬ 
gen  window  purge  to  operate  efficiently.  The  amount  of  ground  propellant 
placed  over  the  ignition  wires  was  reduced  to  minimize  the  effects  of  the 
additional  mass  ejected  during  the  ignition  phase  and  subsequent  clouding  of 
the  windows. 
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Figure  A-2.  Viewing  Section 
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These  changes  significantly  improved  the  initial  photographic  results 
and  made  it  possible  to  observe  the  propellant  surface  during  a  large  portion 
of  the  run.  During  the  early  stages,  a  film  of  opaque  material  was  temporarily 
deposited  or.  the  windows,  preventing  observation  of  the  initial  stages  of  the 
instability  process.  Methods  of  solving  this  problem  are  currently  being 
cons idered . 

A  typical  film  sequence  is  shown  in  figure  A-3  for  the  84 %  AP  propellant. 
The  sides  of  the  test  sample  were  inhibited  before  firing,  and  a  nitrogen 
purge  was  used  in  this  test.  The  photographs  were  taken  from  the  16-mm  color 
film,  and  the  original  film  displays  greater  detail  than  appears  in  figure  A-3. 
Random  frames  were  selected  to  display  various  features  of  the  combustion  pro¬ 
cess.  The  first  frame  shows  the  test  section  before  firing,  and  the  next  frame 
shows  the  ignition  process.  Ignition  is  obtained  by  electrically  heating  the 
nichrome  wire  secured  to  the  propellant  surface,  and  this  appears  as  a  blue 
flash  (i.e.,  bright  flash  near  surface  in  frame  2)  when  the  electrical  contact 
is  broken.  Once  ignition  has  occurred,  mass  is  evolved  from  the  surface,  as 
displayed  by  the  luminous  streaks  near  the  surface  in  frame  3.  At  a  later 
time  the  propellant  started  to  burn  around  the  inhibit.' r,  as  shown  in  frames 
4  and  5 . 


2.0  PHASE  SHIFT  STUDIES 

A  series  of  preliminary  experiments  was  conducted  in  which  the  com¬ 
bustion  process  was  altered  by  phase  shifting  the  pressure  oscillations  with 
electrical  feedback  into  the  combustion  chamber.  When  oscillations  occur,  the 
output  pressure  signal  is  passed  through  a  phase  inverter  and  then  back  into 
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Figure  A-3.  Combustion  Process  in  T-Burner  Viewing  Section 

the  combustion  chamber  through  a  probe  located  inside  the  chamber.  This  tech¬ 
nique  provides  a  possible  method  of  investigating  temperature  oscillations  or 
entropy  waves  in  the  T-burner  and  may  also  moderate  or  eliminate  combustion 
instability  in  motors. 

A  schematic  of  the  feedback  loop  previously  developed  in  connection 
with  flame  acoustic  studies  at  UTC  is  shown  in  figure  A-4.  The  pressure 
transducer  (A)  senses  any  pressure  fluctuations  which  are  then  amplified  (B) 
and  inverted  in  phase  by  180°  (C) .  The  final  signal  is  then  fed  into  the 
c  mbustion  chamber  through  a  probe  (D)  located  above  the  propellant  surface 
and  utilizes  the  case  walls  as  the  ground  electrode. 

The  probe  was  a  tungsten  rod  with  a  ceramic  shield  in  a  connax  fitting. 
It  was  located  about  I  in.  from  the  surface  of  one  of  the  propellant  samples 


A -5 


UTC  2136-FR 
Appendix  A 

in  Che  double-ended  T -burner  and  protruded  into  the  center  of  the  tube. 

The  propellant  used  in  these  preliminary  tests  contained  827,  AP,  2 °L 
CsNO^j  and  16%  UTRE?  hinder.  The  seeding  material,  along  with  the  high  oxidizer 
loading,  was  used  to  inc  ease  the  conductivity  of  the  gases.  The  current,  in 
turn,  was  i elated  to  the  conductivity  and  voltage  across  the  electrodes.  In 
the  present  system  with  500  vdc  across  the  electrode  to  ground,  a  nominal  0.3 
amp  were  obtained  during  the  tests.  As  the  voltage  fluctuated  out-of-phase 
with  the  pressure  due  to  the  feedback  mechanism,  the  current  fluctuated  in-phase 
with  the  pressure,  and  simultaneous  recordings  were  made  of  the  pressure, 
voltage,  and  current  on  the  CEC  recorder. 

These  studies  were  made  at  a  fundamental  oscillating  mode  of  approxi¬ 
mately  350  cps.  Although  it  was  difficult  to  obtain  reproducible  results  from 
these  initial  studies,  several  important  qualitative  observations  were  observed. 
The  results,  based  on  measurements  of  the  chamber  pressure  with  the  probe  in 
the  system  compared  with  and  without  the  500  vdc  and  zero  gain  on  the  feedback 
amplifier,  showed  no  significant  change  in  the  measured  signal.  This  indicates 
that  there  was  no  interference  between  the  electronics  and  the  recorded  pres¬ 
sure.  Small  amplitude  oscillations  appeared  in  the  current  trace  and  were 
in-phase  with  the  pressure  oscillations,  which  is  a  result  of  temperature 
oscillations  that  appear  in  the  T-burner  (i.e.,  '  current  -  -  A  conductivity 
(  '  temperature)^  ). 

Measurements  of  the  chamber  pressure,  made  with  the  probe  in  position, 
were  also  compared  with  and  without  the  feedback  mechanism  in  operation.  The 
results  indicate  that  both  the  absolute  magnitude  in  the  time  varying  pressure 
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Figure  A-4.  Schematic  of  Feedback  Loop 

oscillations  (  A  P)  and  the  sensitivity  of  the  propellant  specified  in  terms 

Af  ,  P 

of  the  acoustic  response  function  (ja/e  =  o  )  were  reduced. 

*o  *  A  P° 
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APPENDIX  B.  EXPERIMENTAL  DATA 


Data  for  the  propellants  tested  or  this  program  are  tabulated  in  this 


appendix.  The  specific  formulations  are  described  in  table  II. 


TABLE  B-I.  UTX-8501  PROPELLANT 

-i 

Density:  0.063  lb/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Freauency 

Growth 
Rate,  */ 

8 

Decay 

Frequency 

Decay 
Rate,  2. 

Re(p/e) 

1 

100  psig 

163 

3.76 

180 

4.48 

2.38 

2 

1 

137 

4.03 

184 

1  .lb 

3.54 

3 

1 

167 

6.18 

193 

4.86 

3.13 

4 

195 

2.35 

205 

11.90 

3.49 

5 

333 

5.30 

320 

1.01 

0.95 

6 

327 

5.88 

296 

0.97 

1.05 

7 

' 

320 

4.53 

291 

2.24 

1.08 

8 

728 

2.20 

667 

1.72 

0.29 

9 

762 

2.16 

762 

2.60 

0.31 

10 

100  i 

psig 

728 

1.81 

635 

2.33 

0.29 

1 

200  psig 

151 

4.72 

151 

2.27 

3.72 

2 

160 

4.24 

144 

3.11 

3.85 

3 

< 

228 

9.42 

242 

3.48 

4.46 

4 

239 

10.20 

266 

2.92 

4.31 

5 

400 

13.47 

219 

7.76 

5.54 

6 

307 

18.81 

285 

2.83 

5.72 

7 

477 

27.40 

785 

5.13 

5.15 

8 

500 

28.55 

741 

7.21 

5.37 

S 

816 

18.77 

579 

11.56 

3.37 

10 

784 

18.42 

615 

10.52 

3.26 

11 

816 

8.98 

563 

15.68 

3.13 

12 

1,904 

13.  ol 

1,818 

6.33 

0.85 

13 

1,812 

13.29 

1,904 

10.04 

1.02 

14 

200  psig 

2,105 

10.34 

1,904 

4.33 

0.57 

1 

500  psig 

140 

2.96 

164 

4.76 

5.20 

2 

186 

3.57 

126 

2.37 

3.95 

3 

177 

1.47 

166 

4.24 

3.51 

4 

226 

9.78 

216 

4.78 

6.75 

5 

220 

10.88 

242 

5.16 

7.32 

6 

260 

12.80 

270 

8.40 

8.33 

7 

762 

7.14 

728 

17.90 

3.52 

8 

667 

5.00 

762 

35.10 

4.17 

9 

o 

o 

isig 

800 

5.30 

762 

34 . 70 

5.40 
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TABLE  B-II .  UTX-8502  PROPELLANT 

Density:  0.063  lb/in.^  Velocity  of  Sound:  34,000  in. /sec 


Run  No . 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  ,.g 

Decay 

Frequency 

Decay 

Rate, 

Re  (\i/c  ) 

1 

200  psig 

7  90 

5.50 

164 

2.23 

3.55 

2 

188 

2.69 

161 

4.72 

3.65 

3 

168 

3.94 

160 

3.20 

3.63 

4 

271 

4.08 

313 

7.4  L 

3.23 

5 

307 

4 . 56 

290 

7.18 

3.30 

6 

242 

4.95 

313 

5.94 

3.28 

7 

667 

a. 49 

635 

3.57 

1.04 

8 

200  psig 

513 

6.77 

769 

1.00 

1.19 

TABLE  B-II1.  UTX-8525  PROPELLANT 

Density:  0.063  lb/in. ^  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  '-g 

Decay 

Frequency 

Decay 
Rate,  ,d 

Re(|i/>-) 

1 

200  psig 

146 

5.18 

131 

4.49 

4.90 

2 

137 

5.00 

139 

5.46 

5.33 

3 

220 

8.59 

223 

9.60 

5.77 

4 

210 

8.43 

218 

8.06 

5.43 

5 

685 

13.47 

647 

3.43 

1.75 

6 

755 

3.21 

631 

13.38 

1.77 

7 

200  psig 

1,620  Could  not  measure  a  growth 

curve  on 

several  runs 
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TABLE  B-IV.  UTX-8526  PROPELLANT 

Density  -  0.063  lb/in^  Velocity  of  Sound  -  36,000  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  'g 

Decay 

Frequency 

Decay 

Rate,:» , 
d 

Re(p/  ■•) 

1 

200  psig 

160 

4.51 

210 

2.35 

2.14 

2 

160 

5.11 

202 

1.98 

2.27 

3 

280 

17.80 

285 

5.37 

4.48 

4 

600 

16.67 

770 

7.35 

2.04 

5 

570 

25.88 

636 

11.74 

3.46 

6 

200  psig 

1500 

Could  not 

measure  a  growth  curve 

on  several  runs 

TABLE  B-V. 

Density  -  0.063  lb/in^ 


UTX-8527  PROPELLANT 

Velocity  of  sound:  36,000  in. sec 


Run  No. 


1 

2 

3 

4 

5 


Mean 

Pressure 


200  psig 


Growth 

Frequency 

Growth 

Rate, 

g 

Decay 

Frequency 

Decay 

Rate, 

Re(  p/f) 

120 

1.23 

165 

1.14 

2.53 

123 

1.41 

170 

1.97 

3.37 

231 

3.59 

275 

3.82 

4.30 

216 

5.76 

263 

2.35 

5.22 

700 

Could  not  measure  a  growth  curve  on 

several  runs 

200  psig 
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TABLE  B-VI.  UTX-8528  PROPELLANT 


3 

Density:  0.063  lb/in.  Velocity  of  Sound:  36,000  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate , 

g 

Decay 

Frequency 

Decay 
Rate,  'd 

Re(  p  /■  ) 

1 

200  psig 

166 

4.95 

146 

5.25 

4.75 

2 

127 

3.64 

143 

5.96 

5.06 

3 

256 

9.72 

245 

2.67 

3.52 

4 

296 

4.28 

256 

6.77 

2.95 

5 

200  psig 

600 

Could  not  measure  a  growth 

curve  on 

several  runs 

TABLE  B-VII.  UTX-8529  PROPELLANT 


Density:  0.063  lb/in.  Velocity  of  Sound:  36,000  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate , 

g 

Decay 

Frequency 

Decay 
Rate,  O !  ^ 

Re(p/c) 

1 

200  psig 

174 

4.65 

176 

5.55 

3.92 

2 

163 

5.20 

161 

4.51 

4.05 

3 

256 

15.29 

285 

4.81 

5.16 

4 

244 

9.79 

284 

6.47 

4.23 

5 

525 

31.32 

730 

5.76 

4.54 

6 

527 

32.70 

784 

6.91 

4.76 

7  200  psig  1,690  Could  not  measure  a  growth  curve  on  several  runs 
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TABLE  B-VIII.  UTX-8531  PROPELLANT 
Density:  0.063  lb/in.  Velocity  of  Sound:  36,000  in. /sec 


Mean 

Growth 

Growth 

Decay 

Decay 

Re(n  /  ) 

Run  No. 

Pressure 

Frequency 

Rate,  ./ 

8 

Frequency  Rate, 

1 

200  psig 

147 

3.66 

167 

5.50 

3.38 

2 

i 

[ 

172 

4.65 

170 

7.00 

3.98 

3 

190 

9.23 

215 

9.16 

5.33 

4 

210 

8.63 

265 

9.10 

4.42 

5 

200 

7.48 

254 

13.88 

5.39 

6 

583 

31.55 

750 

5.85 

3.62 

7 

514 

32.10 

611 

12.73 

4.87 

8 

200  psig 

1,800 

Could  not  mea 

sure  a  growth 

curve  on 

several  runs 

TABLE  B-IX.  UTX-8532  PROPELLANT 


O 

Density:  0.063  lb/in.  Velocity  of  Sound:  36,000  in. /sec 


Mean 

Growth 

Growth 

Decay 

Decay 

Re  (p  /  e) 

Run  No. 

Pressure 

Frequency 

Rate,  -• 

g 

Frequency 

Rate, 

1 

200  psig 

165 

5.83 

135 

6.26 

6.83 

2 

175 

5.69 

136 

5.64 

6.20 

3 

233 

9.83 

232 

7.41 

6.22 

4 

490 

18.31 

755 

15.27 

4.82 

5 

475 

21.51 

690 

4.42 

4.32 

6 

200  psig 

1,500 

3.15 

1,700 

7.11 

0.52 
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TABLE  B-X.  UTX-9164  PROPELLANT 


O 

Density:  0.063  lb/in.  Velocity  of  Sound:  34,000  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 

Ra'.e,  ' 

’  g 

Decay 

Frequency 

Decay 

Rate, 

d 

Re(]i/.  ) 

1 

200  psig 

200 

1.38 

182 

4.26 

2.56 

2 

205 

2.09 

190 

3.02 

2.27 

3 

291 

10.30 

302 

3.41 

3.95 

4 

291 

7.25 

341 

2.60 

2.85 

5 

333 

1.49 

296 

9.10 

2.98 

6 

364 

4.05 

314 

9.57 

3.81 

7 

372 

3.66 

262 

4.07 

2.20 

8 

696 

4.31 

696 

2.43 

0.84 

9 

200  psig 

727 

1.60 

762 

2.38 

0.47 

TABLE  B-XI .  UTX-9165  PROPELLANT 


Q 

Density:  0.063  lb/in.  Velocity  of  Sound:  34,000  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate ,  ' 

g 

Decay 

Frequency 

Decay 
Rate,  ^ 

Re(u  A  ) 

1 

200  p 

sig 

191 

5.11 

170 

3.06 

3.76 

2 

190 

3.66 

174 

3.09 

3.11 

3 

174 

4.93 

174 

2.46 

3.56 

4 

190 

4.46 

170 

3.40 

3.66 

5 

170 

4.62 

178 

1.83 

3.15 

6 

320 

2.43 

250 

3.97 

1.97 

7 

320 

3.30 

222 

2.72 

1.89 

8 

347 

3.06 

235 

3.00 

1.81 

9 

320 

2.23 

242 

5.09 

2.35 

10 

200  psig 

750 

3.24 

750 

2.56 

0.65 
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TABLE  B-XII.  UTX-9167  PROPELLANT 

Density:  0.063  lb/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  it 

E 

Decay 

Frequency 

Decay 

Rate, 

Re(p/ -) 

1 

200  psig 

228 

3.61 

170 

2.67 

1.98 

2 

i. 

216 

3.57 

170 

2.34 

1.90 

3 

222 

4.48 

186 

1.98 

1.94 

4 

223 

3.60 

211 

3.63 

2.10 

5 

333 

7.44 

222 

3.77 

2.46 

6 

320 

4.68 

222 

6.04 

2.70 

7 

800 

3.52 

800 

6.03 

0.75 

8 

200  psig 

800 

3.17 

730 

6.82 

0.84 

TABLE  B-XIII.  UTX-9168  PROPELLANT 


O 

Density:  0.063  lb/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  / 

8 

Decay 

Frequency 

Decay 

Rate, 

Re(p/ c) 

1 

200  psig 

200 

3.36 

163 

5.12 

3.41 

2 

195 

2.96 

157 

4.65 

3.17 

3 

191 

1.69 

145 

5.43 

3.28 

4 

333 

5.38 

250 

7.15 

3.16 

5 

333 

2.58 

210 

6.15 

2.61 

6 

320 

2.65 

265 

6.52 

2.32 

7 

286 

3.24 

308 

6.93 

2.37 

8 

533 

10.95 

730 

7.08 

2.16 

9 

200  psig 

696 

10.20 

800 

3.50 

1.35 
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TABLE  B-XIV.  UTX-9173  PROPELLANT 


Density: 


0.063  lb/ in.3 


Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  ' 

8 

Decay 

Frequency 

Decay 

Rate, 

Re(u/e) 

1 

200  psig 

210 

9.54 

200 

8.60 

5.80 

2 

210 

8.53 

200 

9.71 

5.85 

3 

216 

5.79 

160 

9.80 

5.80 

4 

381 

9.16 

348 

3.98 

2.32 

5 

381 

9.03 

348 

4.16 

2.33 

6 

360 

5.40 

276 

4.58 

2.08 

7 

382 

6.59 

260 

4.91 

2.37 

8 

800 

5.25 

806 

16.56 

1.79 

9 

200  p 

sig 

552 

11.30 

843 

5.88 

1.81 

TABLE  B-XV.  UTX-9368  PROPELLANT 


O 

Density:  0.064  lb/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate , 

g 

Decay 

Frequency 

Decay 
Rate,  a 

P>3(p/e) 

1 

200  psig 

239 

14.78 

267 

2.94 

5.23 

2 

220 

13.38 

267 

2.41 

5.05 

3 

239 

13.21 

276 

3.10 

4.80 

4 

372 

24.90 

762 

8.68 

5.55 

5 

400 

14.75 

727 

7.94 

4.15 

6 

533 

22.70 

800 

8.47 

3.82 

7 

1,010 

6.15 

1,145 

12.40 

1.21 

8 

200  psig 

1,080 

8.32 

1,250 

15.20 

1.43 
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TABLE  B-XV.  UTX-9368  PROPELLANT 


Density:  0.064  lb/in. ^  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  -> 

8 

Decay 

Frequency 

Decay 
Rate,  o:d 

Re  (p  /<: ) 

1 

200  psig 

239 

14.78 

267 

2.94 

5.23 

2 

220 

13.38 

267 

2.41 

5.05 

3 

239 

13.21 

276 

3.10 

4.80 

4 

372 

24.90 

762 

8.68 

5.55 

5 

400 

14.75 

727 

7.94 

4.15 

6 

533 

22.70 

800 

8.47 

3.82 

7 

1,010 

6.15 

1,145 

12.40 

1.21 

8 

200  psig 

1,080 

8.32 

1,250 

15.20 

1.43 

TABLE  B-XVI.  A-13  (NOTS  PROPELLANT) 


Density:  0.057  lb/in. ^  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate,  a 

g 

Decay 

Frequency 

Decay 
Rate,  Ct^ 

Re(|i/e  ) 

1 

200  nsig 

242 

6.4 

157 

10.3 

5.64 

2 

243 

5.4 

182 

13.2 

5.81 

3 

320 

6.5 

314 

13.5 

3.88 

4 

333 

7.2 

286 

8.6 

3.17 

<; 

m/ 

640 

5.2 

572 

8.7 

1.43 

6 

696 

5.5 

552 

10.4 

1.64 

7 

728 

8,9 

762 

12.3 

1.57 

8 

728 

7.1 

572 

11.1 

1.79 

9 

800 

5.2 

800 

17.5 

1.74 

10 

1,430 

11.0 

1,820 

20.6 

1.16 

11 

1,430 

9.5 

1,504 

25.4 

1.44 

<  A 

Lt 

1,540 

14.4 

L  ,  /4u 

13.2 

1.04 

13 

1,670 

10.4 

1,820 

24.6 

1.21 

14 

3,300 

4.3 

3,330 

23.0 

0.51 

15 

J 

' 

3,300 

9.2 

2,850 

37.3 

0.97 

16 

200  p 

-cig 

4,000 

22.3 

4,450 

5.4 

0.42 
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TABLE  B-XVII.  UTX-9379  PROPELLANT 

Density:  0.064  Ib/in."^  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate.  > 

8 

Decay 

Frequency 

Decay 

Rate, 

Re(p/c  ) 

1 

200  psig 

205 

6.45 

154 

2.46 

3.28 

2 

200 

5.79 

167 

1.84 

2.78 

3 

191 

4.66 

170 

1.52 

2.33 

4 

228 

7.62 

165 

4.58 

4.27 

5 

335 

10.22 

281 

4.06 

3.15 

6 

308 

15.3 

262 

2.49 

4.13 

7 

308 

16.16 

291 

3.31 

4.45 

8 

842 

1.16 

800 

3.71 

0.418 

9 

855 

2.58 

855 

3.58 

0.501 

10 

842 

1.18 

842 

2.75 

0.328 

11 

200  psig 

842 

2.67 

762 

2.20 

0.418 

TABLE  B-XVIII.  UTX-9395  PROPELLANT 


Density:  0.064 

lb/ in."* 

Velocity  of  Sound 

36,500 

in. /sec 

i 

Mean 

Growth 

Growth 

Decay 

Decay 

Run  No. 

Pressure 

Frequency 

Rate,  v 

g 

Frequency  Rite,  7, 

Re(p/e  ) 

'  i 

200  psig 

254 

6.93 

210 

10.15 

2.35 

2 

■ 

250 

5.67 

222 

11.63 

2.83 

3 

258 

13.36 

193 

13.02 

3.78 

■  4 

381 

10.27 

286 

7.77 

2.04 

5 

391 

11.71 

391 

8.82 

1.99 

6 

341 

14.22 

334 

9.88 

2.69 

7 

391 

22.85 

308 

5.45 

2.88 

8 

876 

18.65 

800 

44.1 

2.88 

;  ‘  9 

915 

24.5 

855 

60.8 

3.69 

10 

790 

68.65 

627 

10.55 

3.91 

L  11 

842 

29 . 88 

687 

14.8 

2.16 

[  12 

1,880 

16.55 

1,780 

10.00 

0.55 

t  13 

1,780 

14.8 

1,780 

14.95 

0.63 

|  14 

1,780 

23.35 

1,730 

22.00 

0.92 

l  15 

200  psig 

1,780 

11.15 

1,780 

7.7 

0.40 
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TABLE  B-XIX .  UTX-9715  PROPELLANT 

Density:  0.064  lb/in. ^  Velocity  of  Sound:  36  500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate.  og 

Decay 

Frequency 

Decay 
Rate,  a<] 

Re  (n/e) 

1 

200  psig 

228 

10.50 

235 

7.73 

4.56 

2 

f 

266 

6.62 

222 

8.40 

3.63 

3 

242 

7.62 

235 

7.62 

3.85 

4 

276 

6.54 

205 

9.76 

4.14 

5 

235 

7.91 

191 

8.41 

4.50 

6 

336 

7.50 

314 

11.75 

3.47 

7 

381 

7.25 

348 

13.85 

3.42 

8 

348 

9.40 

286 

7.10 

3.01 

9 

76J 

5.24 

696 

3.41 

0.69 

10 

844 

6.27 

667 

3.38 

0.73 

11 

762 

4.20 

842 

8.75 

0.90 

12 

200  psig 

842 

3.37 

842 

2.30 

0.38 

TABLE  B-XX.  UTX-9717  PROPELLANT 

Density:  0.064  lb/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 

Rate ,  cc 
g 

Decay 

Frequency 

Decay 
Rate,  G  d 

ReCu/c) 

1 

200  psig 

239 

5.25 

222 

2.82 

1.67 

2 

203 

9.97 

225 

1.11 

2.55 

3 

232 

5.8 

236 

4.45 

2.12 

4 

372 

3.21 

302 

3.09 

0.91 

5 

364 

4.39 

291 

3.10 

1.10 

6 

364 

6.19 

320 

3.19 

1.29 

7 

381 

3.71 

287 

3.22 

1.58 

8 

854 

4.93 

782 

4.56 

0.559 

9 

728 

4.39 

762 

4.03 

0.545 

10 

842 

1.34 

800 

2.48 

0.505 

11 

200  psig 

842 

4.75 

696 

4.14 

0.559 
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Density:  0.064  Ib/in. 


TABLE  B-XXI.  UTX-9718  PROPELLANT 

1  •  -  3  tt  _  i  . _  c  _ 


Velocity  of  Sound:  36,500  in. /sec. 


Run  No. 


200  psig 


200  psig 


Density:  0.064  Ib/in. 

Mean  Gr 


200  psig 


200  psig 


Growth 

Growth 

Decay 

Decay 

Re(p/e  ) 

Frequency 

Rate,  or 
g 

Frequency 

Rate,  or, 

'  a 

250 

4.62 

296 

3.22 

1.49 

242 

2.38 

286 

1.87 

0.83 

267 

5.0 

281 

4.33 

1.73 

272 

3.23 

239 

3.65 

1.38 

571 

5.79 

421 

5.43 

1.16 

410 

3.75 

364 

3.83 

0.99 

422 

5.08 

422 

7.88 

1.55 

422 

6.94 

410 

3.52 

1.27 

390 

5.91 

258 

4.14 

1.58 

382 

6.58 

291 

5.88 

1.89 

800 

15.65 

667 

6.98 

1.52 

842 

13.53 

941 

7.58 

1.22 

762 

16.0 

616 

3.82 

1.38 

762 

13.15 

640 

5.3 

1.34 

1,730 

17.18 

1,830 

22.45 

1.18 

1,490 

12.1 

1,600 

10.32 

0.73 

1,830 

12.28 

1,940 

11.05 

0.63 

1,780 

9.93 

1,880 

9.75 

0.60 

1,830 

6.9 

1,775 

3.75 

0.25 

TABLE  B-XXII.  UTX- 

9719  PROPELLANT 

b / in. ^ 

Velocity  of 

Sound:  36,500  in. /sec 

Growth 

Growth 

Decay 

Decay 

Re(|i/e) 

Frequency 

Rate,  a  Freouency 

g 

Rate,  ad 

246 

11.47 

281 

7.35 

3.62 

226 

12.27 

235 

9.09 

4.6 

219 

11.74 

250 

6.69 

3.99 

391 

3.12 

281 

3.29 

0.98 

380 

5.83 

286 

3.59 

1.38 

381 

7.08 

281 

5.49 

1.89 

421 

5.38 

421 

7.02 

1.46 

854 

9.83 

705 

5.18 

0.94  • 

865 

14.9 

695 

4.22 

1.18 

842 

14.34 

680 

6.5 

1.33 

1,830 

15.37 

1,560 

'  8.11 

0 .68 

1,939 

21.3 

1,489 

7.65 

0.81 

1,780 

19.55 

1,989 

23.60 

1.14 

1,880 

14.21 

1,680 

4.22 

0.47 
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TABLE  B-XXIII.  UTX-9720  PROPELLANT 

o 

Density:  0.064  Ib/in.  Velocity  of  Sound:  36,500  in. /sec 


Run  No. 

Mean 

Pressure 

Growth 

Frequency 

Growth 
Rate ,  ' 

g 

Decay 

Frequency 

Decay 
Rate ,  «. 
d 

Re(u/(  ) 

1 

200  psig 

254 

2.25 

222 

3.64 

1.35 

2 

254 

3.25 

229 

3.92 

1.58 

3 

211 

3.61 

235 

3.34 

1.66 

4 

411 

7.55 

327 

4.86 

1.78 

5 

372 

7.1 

364 

6.21 

1.9 

6 

390 

5.28 

281 

4.86 

1.63 

7 

372 

5.37 

291 

4.2 

1.52 

8 

762 

16.2 

667 

2.8 

1.4 

9 

800 

14.14 

695 

3.38 

1.18 

10 

800 

12.8 

667 

2.99 

1.08 

11 

890 

10.42 

667 

2.54 

0.82 

12 

1,730 

5.8 

1,640 

4.9 

0.33 

13 

1,775 

10.5 

1,680 

5.65 

0.49 

14 

1,830 

9.35 

1,775 

5.15 

0.15 

15 

200  p 

•sig 

1,775 

6.9 

1,680 

6.3 

0.40 

TABLE  B-XXIV.  UTX-9754  PROPELLANT 

Density:  0.064  lb/in.^  Velocity  of  Sound:  36,500 

Mean  Growth  Growth  Decay  Decay 

Run  No.  Pressure  Frequency  Rate,  o  Frequency  Rate,  / . 

O  ^ 

in. /sec 

Re(p/o) 

1 

200  psig 

188 

4.27 

162 

1.52 

2.9 

2 

154 

5.68 

158 

2.62 

3.93 

3 

156 

4.67 

162 

3.24 

4.52 

4 

145 

6.19 

133 

1.74 

5.04 

5 

286 

3.13 

246 

3.19 

2.16 

6 

239 

6.31 

226 

2.33 

3.33 

7 

250 

5.72 

276 

1.12 

2.43 

8 

...  \ 

271 

9.09 

262 

2.49 

3.9 

9 

200  psig 

577 

4.68 

696 

7.43 

1.68 
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APPENDIX  C 

EQUATIONS  FOR  NONLINEAR  BEHAVIOR 


The  second  order  behavior  of  the  general  analysis  is  given  by 


E(>)  =  Re 


4(14A1-R1) 


t?  o  ^  o  q  ° 

(~)  (— )*  W  +  2  (-i-)  X.  +  Y 

"  1  i  •  i  i  *• 


m  ,  a  o 

'9  ^  i  ...  ^  i 

+  4  (-f)  2  -2  (-i)"  (-M  A 

1  '1  "1  1 


(C-l) 


where 


W  =  (1-R.)[  o2H  fri2)  +  A.  ['-j2?  h$2]  =  (1-R  )A  +  A  R  (C-2) 

1  I  n  Q  I  no  1  1  2 

X  =  (1-R, )  [  o2H  /j?  u  -  ]  +  A  [j2F  /ofio€]  =  (1-R, )C. -  +  A  D  (C-3) 

I  I  1  n  J  o  1  n  Jo  Ill  111 


Yx  =  (1-Rj)  [  ^2HrsA.-:2]0  +  Ai  2]q  =  (l-R^  +  Aj.  B,  (C-4) 


Zx  =  f(l-R1)Q1  +  Ax  BJ 


(C-5) 


However,  by  extrapolating  E(X)  to  E(0)  one  of  the  constants  can  be  eliminated. 
Assuming  =  0  as  X  -»0,  the  resulting  expression  becomes 


E(0) 


(C-6) 


C-l 
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Thus,  can  be  eliminated  from  equation  C-l  yielding  an  equation  which 
co.  tains  three  additional  parameters  (X^,  W^,  and  Z^)  beyond  those  required 
to  characterize  the  linear  region  and  the  steady-state  behavior. 


sc 


The  corresponding  solution  for  ^ /  i  j  the  component  of  which  oscillates 
at  the  second  harmonic,  is 


sc 

^  o2-V:2-V“ 


0  ,  V  2  r 


(Y/V  [W14A2<a2'1>]+  2(61  '  1>  rxi+cil(x2-1)]  +[Y1+Q2(x2_1)3i 
=  -  j 


sc  . 

C2 

+  —^2  j  n(bhA1-R1)+Q1(>2-l)  \  + 


|  _  (u1/e1)(ei°/  x)  I  DmaD-x^-x,) 

2>.1(X1-l)2  M  X2 


(C-7) 


where 


D  =  (J2^‘  l^^l 


m(X)  5  (X1-1)(X2+1-2X1) 


X202-1)  =  iX/2 


The  third  order  behavior  is  given  by  the  equations 
F(>.)  +  A1/>1  -  =  Ax  (A11(X1-1)  +  a12  +  a13  -  a14  -  A17 


\h°\ 

2 

re  ° 

el 

w  t (X1-R)a3+A1R  ] 

r/!i!\2  +  2  A1 

v[(x1-r)c12+a1d1?] 

~1 

/I . 

8  + 

h)  V-i/vJ 

8 

r  m  m 

S+f2 


SC 


t(Xl'Rl)Q2+AlB2] 


C-2 
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Hi  ,  / 6 


o 


+  D1 
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APPENDIX  D 

THEORETICAL  ANALYSIS  OF  SUBSURFACE  REACTIONS  AND  NONLINEAR  EFFECTS 
USING  FINITE  DIFFERENCE  TECHNIQUES 

The  equations  used  to  develop  the  nonlinear  algorithmn  are  essentially 
identical  to  those  being  studied  for  the  linear  behavior.  The  principal 
difference  between  the  two  is  the  manner  in  which  the  subsurface  reaction  term 
is  handled.  In  the  studies  using  the  perturbation  technique,  this  term  is 
neglected;  in  studies  using  the  finite  difference  approach,  this  term  can  be 
considered . 

Using  the  variables  expressed  in  nondimens ional  form 


<?  = 

T/Tqo 

(D-l) 

Cl  = 

v  0 

Ts°/T=o 

(D-2) 

T  = 

t(r  0)2/(> 

(D-3) 

r  ~ 

x  r  In 

O 

(D-4) 

where  r  is  the 
o 

D-l  through  D-4 


steady-state 

become 


combustion  velocity  and  r,  =  k/C 


equa  tions 


’T 


-V 


+  Z„ 
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at  C  =  0 
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Because  of  the  nonlinearity  of  these  expressions,  a  finite  difference  algorithmn 
was  developed  for  digital  computer  solution.  Equation  D-5  in  continuum  form 
replaced  the  finite  difference  approximation,  as  shown  by  equation  D-7,  where 
the  mesh  points  of  the  difference  scheme  are  represented  by  i  from  i  =  1  at  the 
propellant  surface  to  n  at  a  point  nA£  within  the  solid  propellant. 
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Writing  the  boundary  condition  equation  at  t  =  0  in  finite  difference  form 
results  in 


e. 


402  +  301 
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This  equation  may  be  differentiated  with  respect  to  the  time  variable,  t,  to 
give 


(D-9) 


Rearranging 


3  2  A 


-E  /RT  v 
e 


(D-10) 


Thus,  equation  D-7  plus  the  n-1  equations  of  the  form  of  equation  D-4 
provides  a  set  of  n  simultaneous  differential  equations  which  can  be  solved 
with  the  aid  of  a  digital  computer.  Various  expressions  can  be  used  for  P, 
■P/,",  and  >q/jt  and  the  time  differencial  form  of  the  burning  rate  equation 


In  these  studies  a  method  of  solution  has  been  used  which  is  unique  in 
solving  the  transient  diffusion  and  chemical  reaction  equations.  First,  a 
system  of  ordinary  differential  equations  is  generated  from  the  relevant  partial 
differential  equation  by  discretization  in  terns  of  the  space  variable;  then 
the  sequence  of  values  of  the  dependent  variables.  This  leads  to  a  system  of 
ordinary  differential  equations  (in  derivatives  with  respect  to  time).  This 
means  of  solution  of  boundary  value  problems  was  extensively  exploited  for 
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some  time  in  analogue  computer  work  for  solution  of  such  problems  as  transient 
conduction  of  heat  in  solids.  In  analogue  computer  work,  the  limitation  is 
usually  the  capacity  of  the  computer.  The  relatively  recent  development  of 
digital  computer  subroutines  ("package  solvers")  "analogue  discretization" 
methods  have  been  used  as  a  competitive  means  for  solving  certain  types  of 
partial  differential  equations. 

The  specific  advantages  of  this  approach  to  solution  of  boundary  problems 
may  be  listed  as  follows: 

A.  The  programming  time  is  reduced  to  a  minimum  as  most  of  the 
difficult  programming  has  been  done  in  the  development  of  the 
ordinary  differential  equations  subroutine.  Also,  the  computer 
running  time  for  the  differential  equations  approach  seems  to 
offer  an  advantage  compared  to  problems  involving  differencing 
in  both  the  space  and  time  steps. 

B.  Numerical  convergence  problems,  i.e.,  convergence  problems  arising 
in  numerical  solution  of  the  partial  differential  equations,  are 
relegated  to  problems  of  convergence  of  solutions  of  systems  of 
ordinary  differential  equations.  Most  subroutines  for  solution 

of  ordinary  differential  equations  have  some  means  of  assurance  of 
numerical  convergence  of  the  solution, 

C.  Conversion  of  a  boundary  value  problem  into  a  system  of  ordinary 
differential  equations  provides  a  means  of  handling  nonlinear 
difficulties.  This  point  is  particularly  important  because  many 
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of  the  equations  involve  an  exponential  Arrhenius  term  as  well 
as  reactant  product  terms. 

D.  The  capacity  of  the  differenti£ 1  equation  "package  solvers"  is 
sufficiently  large  to  allow  use  of  a  fairly  large  number  of 
mesh  points  in  the  discretization. 

The  study  of  the  acoustic  response  requires  that  the  oscillations  occur 
about  the  steady  state.  Hence,  the  analysis  must  first  predict  the  steady- 
state  conditions  for  a  given  set  of  combustion  parameters  and  then  calculate 
the  transient  behavior.  To  do  this,  two  approaches  are  possible.  One  approach 
is  to  assume  some  initial  temperature  distribution  within  the  solid  and  let 
the  transient  equations  relax  to  the  appropriate  steady  state.  The  pressure 
oscillations  are  then  started  because  the  steady  state  is  reacted. 

The  second  approach  is  to  neglect  ^he  transient  terms  and  calculate  the 
steady-state  conditions  by  a  trial-ana-error  process.  In  this  process,  th^ 
surface  temperature  is  assumed,  and  the  equations  are  solved  using  a  forward 
differencing  procedure  for  the  temperature  deep  within  the  solid.  The  surface 
temperature  is  then  adjusted  until  tin  desired  temperature  in  the  propellant  is 
obtained.  The  resulting  steady-state  p'-ofile  is  then  used  in  the  transient 
mooel,  an>  the  oscillations  are  initi,  ted  immediately. 

Of  these  two  possible  approaches,  the  second  is  preferred  for  economic 
reasons.  The  calculations  using  the  forward  difference  approach  can  be  done 
rapidly  on  the  computer  in  comparison  to  the  transient  calculations.  Hence, 
this  is  the  approach  used  in  these  studies. 
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As  mentioned  in  the  main  body  of  this  report,  even  this  approach  did  not 

prove  to  be  economical.  The  difficulty  arose  when  the  steady-state  profile 

was  inserted  into  the  transient  equations.  Even  though  close  convergence  was 

specified,  the  predicted  steady  state  never  matched  the  exact  steady  state. 

When  such  information  is  inserted  in  the  transient  equations,  further  relaxation 

to  the  true  steady  state  occurs.  This  relaxation  occurs  during  the  initiation 

of  the  oscillations,  and  this  compounding  effect  was  found  to  require  extremely 

2 

small  ratios  of  Ax/Z^  to  maintain  numerical  stability.  This,  in  turn,  results 
in  extremely  long  computer  running  time  (i.e.,  1/2  to  1  hr)  before  the  steady 
oscillatory  behavior  is  achieved.  Many  such  calculations  would  be  required  to 
assess  the  effects  desired;  therefore,  these  studies  were  terminated  because 
of  the  excessive  cost. 
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13.  AIJSTfiACT 
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A  study  has  been  made  to  relc-te  theoretically  predicted  and  experimentally  observed 
pressure-coupled  combustion  instability  behavior  of  solid  propellants.  Using  a 
generalized  theoretical  model  which  incorporates  both  surface  and  gas  phase  combustion 
effects,  the  linear  and  nonlinear  (higher  pressure  amplitude)  region  has  been  character¬ 
ized  by  a  group  of  combustion  parameters.  Concurrent  experimental  studies  have  been- 
directed  toward  investigating  (l)  the  effects  of  coating  the  solid  oxidizer  crystals, 

(2)  the  chemical  and  physical  nature  of  the  propellant  binder,  catalysts,  and  cross- 
linking  agents  on  the  acoustic  admittance,  (3)  nonlinear  response  properties,  and 
00  propellant  burning  rate.  The  data  show  that  small  changes  in  propellant  formulatic; 
can  have  a  significant  effect  on  the  derived  parameters  cha ? acteriz ing  the  surface  and, 
j therefore,  on  relative  contributions  of  these  reactions  to  the  acoustical  response 
properties. 

Higher  order  effects  show  that,  depending  on  the  magnitude  of  the  constants,  both 
positive  and  negative  values  for  the  shift  in  the  time  average  burning  rate  with  acous. 
j pressure  can  be  obtained.  Additional  studies  show  that  the  acoustic  response  function 
decreases  with  increasing  acoustic  pressure  for  reasonable  second  and  third  order 
coefficients.  Fence,  tine  combustion  process  itself  can  15-.' t  the  acoustic  pressure 
amplitude  in  the  T- turner.  On  the  oasis  of  these  calculations,  limiting  amplitudes  of 
,0.1  to  0.5  can  be  expected  without  consideration  of  acoustic  losses.  The  combined 
•linear  and  nonlinear  results  indicate  that  a,  homogeneous  nonreacting  solid  character 5 
sthe  transient  behavior  of  the  solid  phase  with  sufficient  accuracy. 
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